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ON LIPSCHITZ SOLUTIONS FOR SOME 
FORWARD-BACKWARD PARABOLIC EQUATIONS 

SEONGHAK KIM AND BAISHENG YAN 


Abstract. We investigate the existence and properties of Lipschitz so¬ 
lutions for some forward-backward parabolic equations in all dimen¬ 
sions. Our main approach to existence is motivated by reformulating 
such equations into partial differential inclusions and relies on a Baire’s 
category method. In this way, the existence of infinitely many Lipschitz 
solutions to certain initial-boundary value problem of those equations 
is guaranteed under a pivotal density condition. Finally, we study two 
important cases of forward-backward anisotropic diffusion in which the 
density condition can be realized and therefore the existence results fol¬ 
low together with micro-oscillatory behavior of solutions. The first case 
is a generalization of the Perona-Malik model in image processing and 
the other that of Hollig’s model related to the Glausius-Duhem inequal¬ 
ity in the second law of thermodynamics. 


1. Introduction 

The evolution process of many quantities in applications can be modeled 
by a diffusion partial differential equation of the form 

(1.1) ut = div{A{Du)) innx(0,T), 

where 11 C M"' is a bounded domain, T > 0 is any fixed number, and 
u = u{x, t) is the density of the quantity at position x and time t, with 
Du = (uxi, • • • j Uxn) Ut denoting its spatial gradient and rate of change, 
respectively. The vector function A: MA —)• MA here represents the diffusion 
flux of the evolution process. The usual heat equation corresponds to the 
case of isotropic diffusion given by the Fourier law: A{p) = kp {p € M”), 
where fe > 0 is the diffusion constant. 

For standard diffusion equations, the flux A{p) is assumed to be monotone; 
namely, 

{A{p) - A(q)) • (p - g) > 0 (p, q€ M”). 

In this case, equation 0 is parabolic and can be studied by the standard 
methods of parabolic equations and monotone operators P EH Eo]. In 
particular, when A(p) is given by a smooth convex function W{p) through 
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A{p) = DpW{p) {p E W^), 113 can be viewed and thus studied as a certain 
gradient flow generated by the energy functional 


I{u) = / W{Du{x))dx. 

Jn 

However, for certain applications of the evolution process to some im¬ 
portant physical problems, the underlying diffusion fluxes A(p) may not be 
monotone, yielding non-parabolic equations (1.1). In this paper, we study 
the diffusion equation (1.1) for some non-monotone diffusion fluxes A{p). 
We focus on the initial-boundary value problem 


( 1 . 2 ) 


ut = div(H(Ilu)) in Qt, 

< A{Du) • n = 0 on dO, x (0,T), 
u = uq on X {t = 0}, 


where Ut = H x (0,r), n is the outer unit normal on 911, uq = uo(x) is a 
given initial datum, and the flux A(p) is of the form 

(1.3) A(p) = f{\p\‘^)p {p E M"), 

given by a function /: [0, oo) —)> M with profile a{s) = sf{s‘^) having one of 
the graphs in Figures and below. (Precise structural assumptions on 
it(s) = s/(s^) will be given in Section]^) 



Figure 1. Case I: Perona-Malik type of profiles a{s). 

The two cases in Figures[^[^and[^correspond to the applications in image 
processing proposed by Perona and Malik ^5] and in phase transitions of 
thermodynamics studied by Hollig [21]. For these diffusions, we have 

it'(\/s) = f{s) + 2s/'(s) < 0 for some range of s > 0. 

In these cases, the diffusion is anisotropic since the diffusion matrix {Ap.), 
where 


^p, (p) =/(IPp)^ii + 2/'(|pP)piPj (i, j = 1,2, • • • ,n). 
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Figure 2. Case II: Hollig type of profiles (t(s) with 
(t(s2) > 0. 



Figure 3. Case II: Hollig type of profiles (t(s) with 
cr(s2) = 0. 


has the eigenvalues /(|pp) (of multiplicity n — 1) and /(|pp) + 2|pp/'(|pp); 
hence the diffusion coefficients could also be negative, and problem (1.2) 
becomes forward-backward parabolic. Moreover, setting 


r\p\ 

W{p)= / a{r)dr, 

Jo 


I{u) = / W{Du) dx, 
Jn 
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the initial-boundary value problem (1.2) becomes a L^-gradient flow of the 
energy functional I{u); however, I{u) is non-convex. Consequently, neither 
the standard methods of parabolic equations and monotone operators nor 


the non-linear semigroup theory can be applied to study (1.2) 


Forward-backward parabolic equations have found many important ap¬ 
plications in the mathematical modeling of some physical problems, but 
mathematically, due to the backward parabolicity, the initial-boundary value 


problem (1.2) of such equations is highly ill-posed, and in many cases even 


the notion and existence of reasonable solutions remain largely unsettled. 

In |35j . the original Perona-Malik model (1.2) in image processing used 
the profile 


(1.4) 


a 


(s) = se or ct(s) = 


l + s^sl 


for denoising and edge enhancement of a computer vision; in this case, we 


call (1.1) the Perona-Malik equation. In this model, u{x,t) represents an 
improved version of the initial gray level uq{x) of a noisy picture. The 
anisotropic diffusion div{A(Du)) is forward parabolic in the subcritical re¬ 
gion where \Du\ < sq and backward parabolic in the supercritical region 
where \Du\ > sq- The expectation of the model is that disturbances with 
small gradient in the subcritical region will be smoothed out by the forward 
parabolic diffusion, while sharp edges corresponding to large gradient in the 
supercritical region will be enhanced by the backward parabolic equation. 
Such expected phenomenology has been implemented and observed in some 
numerical experiments na [35], showing the stability and effectiveness of 
the model. Mathematically, there have been extensive works on the Perona- 
Malik type of equations with profiles (j{s) as in Figure however, most of 
these works have focused on the analysis of numerical or approximate solu¬ 
tions by different methods. For example, in dimension n = 1, the singular 
perturbation and its F-limit related to the Perona-Malik type of equations 
were investigated in 13 El- In |2n) . a mild regularization of the Perona- 
Malik equation with a viscous term was used to extract a unique approxi¬ 
mate Young measure solution in any dimension. The works Pl39| studied 
Young measure solutions for the Perona-Malik equation in dimensions n = 1 
and n = 2. Also the works [I5lll6| focused on the numerical schemes for 
the Perona-Malik model. Recently, classical solutions for the Perona-Malik 
equation were studied in [El [El, where the existence of solutions was proved 
for certain initial data uq that can be transcritical in the sense that the sets 
{\Duo{x)\ > So} and {|L>tto(x)| < sq} are both nonempty (sq is the number 
in Figure [^; however, their initial data cannot be arbitrarily prescribed. 

For the Hollig model, the one-dimensional forward-backward parabolic 
problem ( |1.2[) under the piecewise linearity assumption on the profile a{s) 
(see Figure was studied in mi, motivated by the Clausius-Duhem in¬ 
equality in the second law of thermodynamics in continuum mechanics. For 
such a special profile, it was proved that there exist infinitely many L^-weak 
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solutions to in dimension n = 1. The piecewise linearity of a{s) in 

n = 1 was much relaxed later to include a more general class of profiles 
it(s) (as in Figure]^ in the work of Zhang [H], using an entirely different 
approach from Hollig’s. 



The question concerning the existence of exact weak solutions to problem 


(1.2) of the Perona-Malik and Hollig types had remained open until Zhang 


[la HU first established that the one-dimensional problem ( |1.2[ ) of each 
type has infinitely many Lipschitz solutions for any suitably given smooth 
initial data uq. Zhang’s pivotal idea was to reformulate equation (1.1) in 
n = 1 into a 2 X 2 non-homogeneous partial differential inclusion and then 
to prove the existence by using a modified method of convex integration, 
following the ideas of |28[ [32] • The study of general partial differential 
inclusions has stemmed from the successful understanding of homogeneous 
differential inclusions of the form Du{x) G K, first encountered in the study 
of crystal microstructures la El [32]; recently, the methods of differential 
inclusions have been successfully applied to other important problems; see, 
e.g., iini[iaEiiE3iEsi[ii]. 

In general, a function u G is called a Lipschitz solution to 


problem (1.2) provided that equality 


(1.5) / {u{x, s)({x, s) — uo{x)({x,0))dx = / / {uQt — ^{Du) ■ DQdxdt 

Jn Jo Jn 

holds for each G 6*°° (fir) and each s G [0, T]. Let (" = 1; then it is im¬ 
mediate from the definition that any Lipschitz solution u to (1.2) conserves 
the total mass over time: 


/ u{x,t)dx = / uo{x)dx VtG[0,T]. 
In Jn 
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In [26], we extended Zhang’s method of partial differential inclusions 
to the Perona-Malik type of equations in all dimensions n for balls = 
Br{ 0) := {x G M"’ I |x| < -R} and non-constant radially symmetric smooth 
initial data uq. In this case, the n-dimensional equation for radial solutions 
can still be reformulated into a 2 x 2 non-homogeneous partial differential 
inclusion, and the existence of infinitely many radial Lipschitz solutions to 
(1.2) is established. However, for general domains and initial data, the 


n-dimensional problem (1.2) can only be recast as a (1 -|- n) x (n -|- 1) non- 
homogeneous partial differential inclusion that has some uncontrollable gra¬ 
dient components, making the construction of Lipschitz solutions for this 
inclusion impractical. In our recent work 1211, we overcame this difficulty 
by developing a suitably modified density method, still motivated by the 
method of differential inclusions but based on a Baire’s category argument. 
In this work, it was proved that for all bounded convex domains H C M” with 
cIH G and initial data uq G with Duq-h = 0 on (9H, there exist 

infinitely many Lipschitz solutions to (1.2) for the exact Perona-Malik dif¬ 
fusion flux A{p) = {p G M"). However, the proof heavily relies on the 

explicit formula of the rank-one convex hull of a certain matrix set defined 
by this special flux; such an explicit formula for the general Perona-Malik 
type of flux A{p) (with profile a{s) as in Figure is unattainable. 

The main purpose of this paper is to generalize the method of [27| to 
problem 0 in all dimensions for general diffusion profiles of the Perona- 
Malik and Hollig types. To state our main theorems, we make the following 
assumptions on the domain H and initial datum uq: 


( 1 . 6 ) 


H C M"" is a bounded domain with dfi of 

Uq G is non-constant with Duq ■ = 0, 


where a G (0,1) is a given number. Without specifying the precise con¬ 
ditions on the profiles (t(s) for Cases I and II, we first state our main 
existence theorems as below. 


Theorem 1.1. Assume condition (1.6) is fulfilled. Let A{p) be given by 


(1.3) with the profile cr{s). Then problem (1.2) has infinitely many Lipschitz 
solutions in the following two cases. 

Case I: a{s) is of the Perona-Malik type as in Figur^^ and in addition 
Q is convex. 

Case II: cr(s) is of the Hollig type as in Figur^^ and in addition uq 
satisfies |Zlrio(2:o)| G (s^,^^) at some xq G H. 

The precise assumptions on the profiles a{.s) and detailed statements of 


the theorems for both cases will be given later (see Theorems 3.1 and 3.5) 


along with some discussions about a certain implication of our result on the 
Perona-Malik model in image processing and the breakdown of uniqueness 
for the Hollig type of equations. 
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In certain cases, among infinitely many Lipschitz solutions, problem (1.2) 
admits a unique classical solution; we have the following result (see also 


Theorem 1.2. Assume condition (m is fulfilled and 11 is convex. Let A{p) 
he given by (1.3) with the profile (t{s). Assume ||Z1uo||loo(q) < sq if (x{s) is 


of the Perona-Malik type or ||I1mo||l°°(o) < ■si if <^{s) is of the Hollig type. 
Then problem (1.2) has a unique solution u G satisfying 

This theorem will be also separated into Cases I and II in Section in 
order to be compared with our detailed main results. Theorems 3.1 and |3.5[ 
Finally, we have the following coexistence result. 

Theorem 1.3. Assume that Ll = is a ball in M” and that uq is 

radial and satisfies condition (1.6). Then there are infinitely many radial 


and non-radial Lipschitz solutions to (1.2). 


The rest of the paper is organized as follows. Section begins with 


a motivational approach to problem (1.2) as a non-homogeneous partial 


differential inclusion and its limitation. Then the general existence theorem, 
Theorem 2.2 is formulated under a key density hypothesis, and some general 
properties of Lipschitz solutions to (1.2) are also provided. In section]^ 


precise structural conditions on the profiles a{s) of the Perona-Malik and 
Hollig types are specified, and the detailed statements for the main result, 
Theorem [13 are introduced as Theorems |3.1| and |3.5| that are followed by 
some remarks and related results. Section prepares some useful results 
that may be of independent interest. As the core analysis of the paper, 
the geometry of related matrix sets is scrutinized in Section leading to 
the relaxation result on a homogeneous differential inclusion, Theorem |5.8[ 
Section is devoted to the construction of suitable boundary functions and 


admissible sets for Cases I and II in the contexts of Theorems 13.11 and 3.5 


respectively. Then the pivotal density hypothesis in each case is realized in 
Section completing the proofs of Theorems ^ and ^ Lastly, Section 
deals with the proof of the coexistence result, Theorem |1.3[ 

Throughout the paper, we use the boldface letters for Cases I and II for 
clear distinction. The paper follows a parallel exposition to deal with both 
cases simultaneously. But for a better readability, we strongly recommend 
for a reader to follow each case one at a time. 


2. Existence by a general density approach 

Here and below, we assume without loss of generality that the initial 
datum Mo £ in problem (1.2) satisfies 


uq{x) dx = 0, 


(2.1) 
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since otherwise we may solve (1.2) with initial datum uq = uq — uq, where 
'“0 ~ PI /n ^0 dx. 

2.1. Motivation: non-homogeneous partial differential inclusion. 


To motivate our approach, let us reformulate problem (1.2) into a non- 
homogeneous partial differential inclusion. 

First, assume that we have a function = {u*,v*) G 
satisfying 

u*{x,0) = uo{x), X G n, 

(2.2) < diYV*{x,t) = u*{x,t), a.e. {x,t) ^ ^t, 

• n|ao = 0, t€[0,T], 

which will be called a boundary function for the initial datum uq. 

Let A G be the diffusion flux. For each I G M, let K{1) be a 

subset of the matrix space defined by 


(2.3) K{1) = 


P 

B 


A{p) 


p e 


c G 


B G 


tr B = I 


Suppose that a function w = {u, v) G solves the Dirichlet 

problem of non-homogeneous partial differential inclusion 


(2.4) 


Vw{x,t) ^ K{u{x,t)), a.e. {x,t)^QT, 
w{x,t) = ^(x,t), {x,t) £ d^T, 

where Vw is the space-time Jacobian matrix of w: 

^ Du ut' 

Dv vt 


Vtc = 


Then it can be easily shown that u is a Lipschitz solution to (1.2). 


Although solving problem (2.4) is sufficient to obtain a Lipschitz solution 


to (1.2), all known g enera l existence results [laiM are not applicable to 
solve (2.4). If n = 1, (2.4) implies Vx = u and thus v G if u G 


so Zhang [431144] was able to solve (2.4) for the Perona-Malik and Hollig 


types in dimension n = 1 by using a modified convex integration method. 
However, for n > 2, since it is impossible to bound \\Dv\\loo(^q^^ in terms of 
divu, the function v may not be in K”") even when u G 

therefore, solving (2.4) in M^+”') is impractical in dimensions n > 


2. To overcome this difficulty, we make the following key observation. 


Proposition 2.1. Suppose u G is such that tt(x,0) = uo{x), 

there exists a vector function v G VF^’^((0, T); L^(H; M"’)) with weak time- 
derivative Vt satisfying 

(2.5) Vt = A{Du) a.e. in^T, 

and for each ( G and each t G [0, T], 

(2.6) / v(x,t) ■ D({x,t) dx = — / u{x,t)C{x,t) dx. 

Jo. J o 
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Then u is a Lipschitz solution to (1.2). 


Proof. To verify (1.5), given any C £ let 

g{t) = / u{x,t)({x,t)dx, h{t) = / u{x,t)<ft{x,f)dx (tG[0,T]). 

J ^ Jo. 

for each G Cf°{0,T), 

Tptv ■ D(dxdt, / fjhdt = — 'ipv ■ DQdxdt. 

Jo Jo Jn 


Then by 

ri;,gdt = - r j 

Jo Jo Jn 

Since v G 1T^’^((0, T); L^(ri; M”)) and vt = A{Du) a.e. in Qt, we have 


10 Jn 


{fjDQt • vdxdt = — f f A{Du) ■ ifDf dxdt. 

Jo Jn 


As {f)DQ)t = i^tDC, + combining the previous equations, we obtain 

f dt = [ 'tf(—h+ [ A{Du) ■ D( dx'\ dt, 

Jo Jo \ Jn J 

which proves that g is weakly differentiable in (0, T) with its weak derivative 

g'{t) = h{t)— / A{Du{x,t)) ■ Df{x,t) dx a.e. tG(0, T). 

Jn 

From this, upon integration, (1.5) follows for each s G [0,T]. □ 

2.2. Admissible set and the density ap proach. Let = {u*,v*) be 
any boundary function for uq defined by (2.2) above. Denote by 

M"") the usual Dirichlet classes with boundary traces n*, v*, re¬ 
spectively. 

We say that U C is an admissible set provided that it is 

nonempty and bounded in and that for each u £U, there exists 

a vector function v G M"') satisfying 

divu = u a.e. in D-r, ||ut||£,oo(Qy) < i?, 

where i? > 0 is any fixed number. If U is an admissible set, for each e > 0, let 
Ue be the set of all u G such that there exists a function v G 1^”) 

satisfying 

divu = u a.e. in < R, 


/ \vt{x, t) — A{Du{x, t))| dxdt < elD^I. 

J 

The following general existence theorem relies on a pivotal density hy¬ 
pothesis of Uf: in U. 


Theorem 2.2. Let lA C {TIt) be an admissible set satisfying the den¬ 

sity property: 

(2.7) Ue is dense in lA under the L°^-norm for each e > 0. 
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Then, given any (/? £ U , for each 6 > 0, there exists a Lipschitz solution 
u G to (1.2) satisfying ||m — (p\\ L°°{nT) < Furthermore, if U 

contains a function which is not a Lipschitz solution to (1.2), then (1.2) 
itself admits infinitely many Lipschitz solutions. 

Proof. For clarity, we divide the proof into several steps. 

1. Let X be the closure of Id in the metric space L°°{Ll'j')- Then (T, L°°) 

is a nonempty complete metric space. By assumption, each is dense in 
X. Moreover, since Id is bounded in we have X C 

2. Let y = For h > 0, define Th'. X ^ y as follows. Given 

any u G T, write u = u* + w with w G and define 

Thiu) = Du* + D{ph * w), 

where ph{z) = h~^p{z/h), with z = {x,t) and iV = n + 1, is the standard 
/i-mollifier in , and Ph*w is the usual convolution in with w extended 
to be zero outside Then, for each h > 0, the map : {X, L°°) —>■ (T, L^) 
is continuous, and for each u G X, 

lim \\Th{u) - DuWl^^ ) = lim \\ph * Dw - Dw\\l^q ^ = 0. 

h^0+ ' ' h^0+ 

Therefore, the spatial gradient operator D: X ^ y is the pointwise limit of 
a sequence of continuous maps T^: X ^ y-, hence D: X ^ y is a Baire-one 
map. By Baire’s category theorem (e.g., [3 Theorem 10.13]), there exists a 
residual set Q <Z X such that the operator D is continuous at each point of 
Q. Since T \ ^ is of the first category, the set Q is dense in X. Therefore, 
given any </? G T, for each 6 > 0, there exists a function u G Q such that 
11 ^ ~ pWl^XIt) < 


3. We now prove that each ir G ^ is a Lipschitz solution to (1.2). Let u G Q 


be given. By the density of Id,, in {X, L°°) for each e > 0, for every j G N, 


there exists a function Uj G Idxjj such that 


n,- — tt 


< 1/j- Since the 


operator D: {X,L°°) -G {y,L^) is continuous at u, we have Duj -G Du in 


w^r{nT-, 


L^{LIt','^^). Furthermore, from (2.2) and the definition of ^i/j, there exists 
a function Vj G 
t G [0, T], 


such that for each ( G and each 


( 2 . 8 ) 


Vj{x,t) ■ D({x,t) dx = — / Uj{x,t)Q{x,t) dx, 

! Jn 


\\ivj)t\\L°°{nT) — F, 


\{vj)t — A{Duj) \ dxdt < -|nr|- 


Since Vj{x,0) = n*(x,0) G lF^’°°(n;M"') and \\{vj)t\\L°°{nT) — F, it fol¬ 
lows that both sequences {vj} and {(iij)*} are bounded in L‘^{Qt',^"') ~ 
L^((0,r);L^(n;M"')). So we may assume 

Vj V and {vj)t vt in L^((0, T); M"')) 

for some v G VF^’^((0, T); M"')), where ^ denotes the weak conver¬ 

gence. Upon taking the limit as j —)• oo in (2.8), since i; G C([0, Tj; M"-)) 
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and A € we obtain 

/ v{x,t) ■ D({x,t) dx = — / u{x,t)C{x,t) dx (tG[0,T]), 
JQ. Jq 


vt{x,t) = A{Du{x,t)) CL.e. {x,t) £ ^^r- 


Consequently, by Proposition |2.1[ n is a Lipschitz solution to (1.2). 

4. Finally, assume U contains a function which is not a Lipschitz solution 
to (1.2); hence Q ^U. Then Q cannot be a finite set, since otherwise the 
L°°-closure X = Q = U would be a hnite set, making U = Q. Therefore, 

■is 

□ 


in this case, (1.2) admits infinitely many Lipschitz solutions. The proof is 
complete. 


The following result implies that the density approach streamlined in 
Theorem 2.2 can be useful only when problem (1.2) is non-parabolic, that 
is, when A{p) is non-monotone. 

Proposition 2.3. If A: 

one Lipschitz solution. 

rl.oo 


is monotone, then (1.2) can have at most 


Proof. Let u, u £ be any two Lipschitz solutions to (1.2). Use the 

identity (1.5) for both u and u, subtract the two equations, take C = u — u 
(by approximations), and then apply the monotonicity of A{p): 

{A{Du) — A{Du)) ■ {Du — Du) > 0 

to obtain 

/ {u{x, s) — u{x, s))‘^dx < / / {u — u){u — u)tdxdt 

.Jo, Jo Jn 

= - I {u{x, s) — u{x, s))"^dx (0 < s < T). 

2 Jn 

Thus u = uva. Qt- D 


As a general property for Lipschitz solutions to (1.2) for all continuous 
fluxes A{p) satisfying a positivity condition below, we prove the following 
result; clearly, the positivity condition is satisfied by the fluxes A{p) given by 
(0 with the profiles a{s) of the Perona-Malik and Hollig types illustrated 
in Figures [T} [2] and Note that the positivity condition is consistent with 
the Clausius-Duhem inequality in the second law of thermodynamics (e.g., 
dal p. 79] and [lOl p. 116]). 


Proposition 2.4. LetuQ G VF^’°°(f2) and A: 
condition: 

A{p) ■ p>0 (p G I 


satisfy the positivity 


Then any Lipschitz solution u to (1.2) satisfies 

(2.9) minuo If u ^ maxtto in Dt- 

n o 
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Proof. Let u G 
all C e 


{^t) be any Lipschitz solution to (1.2). By (1.5), for 



ut{x, t)C{x, t)dxdt 



A{Du) ■ D(dxdt; 


hence by approximation, this equality holds for all C £ Taking 

C{x,t) = (f{x,t)ilj{t) with arbitrary f G W^'^{PtT) and if G W^’°°{0,T), we 
deduce that 


ut{x,t)(t>{x,t) dx = — / A{Du{x,t)) ■ D4>{x,t) dx 

! Jn 


for a.e. t G (0,r) and all (f G Now taking cf = ^ with 

A; = 1, 2, • • •, we have for a.e. t G (0, T), 


fL 

dt 


I u^^di^ = 2k I Ut4>dx = —2k / A{Du) ■ Dfdx 

Jq ^ Jq Jq 

= -2k{2k - 1) / u^^-‘^A(Du) ■ Dudx < 0. 

Jn 


From this we deduce that L^^-norm ofu{-,t) is non-increasing on t G [0,T]; 
in particular, 


^)llL2fe(n) < ||^o||L2fc(f^) VtG[0,T], k — l,2,---. 
Letting k —>• oo, we obtain ||u(-,< ||uo||loo(q); hence 
(2.10) = ||tto||L°°(Q)- 


Now let mi = min^ uq and m 2 = niax^ uq- We show mi < u(x, t) < m 2 for 
all (x, t) G to complete the proof. We proceed with three cases. 

(a): m 2 > 0 and \mi\ < m 2 . In this case, ||uo||loo(q) = m 2 ] so by 


( 2 . 10 ) 


u{x,t) < ||n||i,oo(Qj,) = ||uo||L°°(n) = "i 2 - 


To obtain the lower bound, let uq = —uq + m 2 + mi and u = —u + m 2 + mi. 
Then u is a Lipschitz solution to (1.2) with new flux function A{p) = —A(—p) 
and initial data uq. Since mi < uo{x) < m 2 , as above, we have u{x, t) < m 2 ] 
hence u{x,t) > mi for all {x,t) G Ht- 

(b) : m 2 > 0 and mi < —m 2 . Let uq = —uq and u = —u. Then u is a 
Lipschitz solution to (1.2) with new flux function A{p) = —A{—p) and initial 
data Uq. Since —m 2 < uo{x) < —mi for all x G 11 and —mi > 0, | — m 2 \ = 
m 2 < —mi, it follows from Case (a) that —m 2 < u{x,t) < —mi and hence 
mi < u{x,t) < m 2 for all (x, f) G Ht- 

(c) : m 2 < 0. In this case mi < 0. If mi = 0 then m 2 = 0 and hence 
Uq = 0; so, by (2.10), u = 0. Now assume mi < 0. Let again as in Case 
(b) Uq = —Uq and u = —u. Since —m 2 < uq{x) < —mi for all x G 11 and 
—mi > 0, I — m 2 1 = —m 2 < —mi, it follows again from Case (a) that 
—m 2 < u{x,t) < —mi and hence mi < u{x,t) < m 2 for all (x,l) G Ht- TI 
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The rest of the paper is devoted to the construction of suitable boundary 
functions = {u*,v*) and admissible sets U C Wlf°{QT) fulhlling the 
density property (2.7) for Cases I and II. 


3. Structural conditions on the profiles: 

DETAILED STATEMENTS OF MAIN THEOREMS 


In this section, we assume the domain and initial datum uq satisfy (1.6) 
We consider the diffusion fluxes A{p) given by (IL 


and present the detailed 
statements of our main theorems by specifying the structural conditions on 
the profiles a{s) = sf{s‘^) illustrated in Figures and 


3.1. Case I: Perona-Malik type of equations. In this case, we assume 
the following structural condition on the profile cr{s). 

Hypothesis (PM): (See Figuresand [^) 

(i) There exists a number sq > 0 such that 

/ G c°([o, oo)) n c3([o, s§)) n oo). 

(ii) cr'(s) > 0 Vs G [0,So), cr'is) < 0 Vs G (so,oo), and 

lim cr(s) = 0. 

S^OO 

In this case, for each r G (0,cr(so)), let s_(r) G (0, sq) and s+(r) G (so,oo) 
denote the unique numbers with r = cr(s±{r)). Then by (ii). 


(3.1) 


lim s_(r) = 0, lim S-ufr) = oo. 
)—>-0+ r—>•0+ 


Note that both profiles in (1.4) for the Perona-Malik model 
Hypothesis (PM). 

The following is the first main result of this paper in detail. 


satisfy 


Theorem 3.1. Let Ll be convex and Mq = ||Zlno||ioo(Q). Then for each 
r G (0, cj(Mo)), there exists a number I = Ir G (0, r) such that for all r G {I, r) 
and all but at most countably many r G (0,f), there exist two disjoint open 
sets C LIt with U = {LItI and infinitely many Lipschitz 

solutions u to (1.2) satisfying 

u G C'^“'““’^“'“"/^(Hy), ut = div{A{Du)) pointwise in 


\Du{x,t)\ < s_(r) V(x, f) G Hq C 

IS"! -|- \L\ = |H^|, \L\ > 0, 

where 

Hg = {(x, 0) I X G |T>uo(x)| < s_(f)}, 

S = {(x,f) G I s_(f) < \Du{x,f)\ < 'S-(f)}, 

and 

L = {(x,t) G I s+(r) < \Du{x,t)\ < s+(r)}. 
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Remark 3.2. As we will see later, once the numbers f, f are chosen as 
above, the corresponding Lipschitz solutions in the theorem are all identi¬ 
cally equal to a single function u* in which is the classical solution to 
some modified uniformly parabolic Neumann problem. Here the function u* 
depends on f but not on f, whereas do on f: 

= {(®) t) e Qt\ \Du*{x, t)| < s_(f)}, 

= {(x, t) G I \Du*{x, t)| > 

The choice of f is made to guarantee that the interface Qt \ ^r) 

has (n -|- l)-dimensional measure zero; this will be crucial for the proof of 
the theorem later. Choosing r = r may not be safe in this regard, since we 
have not enough information on the function u* to be sure that the interface 
measure |{|T)tt*| = s_(f)}| = 0. This forces us to sacrifice the benefit of the 
choice f = r that would separate the space-time domain into two disjoint 
parts where the Lipschitz solutions are ( 72 +«,i+a /2 nowhere 

in the other. 


Remark 3.3. By (3.1), if 0 < r <C a{Mo), the corresponding Lipschitz 
solutions u have large and small gradient regimes L and U S in Qt up 
to measure zero, representing the almost constant and sharp edge parts of u 
in Qt, respectively. Although there is a fine mixture of the disjoint regimes 
L, S C due to a micro-structured ramping with alternate gradients of 
finite size, such properties together with ( |2.9[ ) for solutions u are somehow 
reflected in numerical simulations; see Figure taken from Perona and Ma¬ 
lik [3^. On the other hand, it has been observed in 011] that as the limits 
of solutions to a class of regularized equations, infinitely many different evo¬ 
lutions may arise under the same initial datum uq. Our non-uniqueness 
result seems to reflect this pathological behavior of forward-backward prob¬ 


lem (1.2) 


We now restate Theorem |1.2| only for Case I. The proof of this appears 
in Subsection! 


The orem 3.4. Let H and Mq be as in Theorem 3.1. If Mq < sq, then 
1.2) has a unique solution u G satisfying \\Du\\icx>(^q.^'j = 

(n)- 


'm 




For the given initial datum uq with Mq < sq, by Theorem 3.1, there are 


infinitely many Lipschitz solutions to problem (1.2). On the other hand, by 


the theorem here, we have a special Lipschitz solution u to (1.2), which is 


also classical. It may be interesting to observe that none of the Lipschitz 


solutions from Theorem 3.1 coincide with this special solution u. 


3.2. Case II: Hdllig type of equations. In this case, we impose the 
following structural condition on the profile ct(s). 

Hypothesis (H): (See Figures]^ and [^) 
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(a) 


z 



0.0 


(b) 


Figure 5. Scale-space using anisotropic diffusion profile 
~ i_|_^ 2 /g 2 ■ Three dimensional plot of the brightness of 

Fig. 12 in [35]. (a) Original image, (b) after smoothing with 
anisotropic diffusion. 


(i) There exist two numbers S 2 > si > 0 such that 


/ G C°([0, oo)) n C'^+"([0, sf) U (si oo)). 


(ii) cr'(s) > 0 Vs G [0,si) U (s 2 ,oo), (t(si) > ( 7 ( 82 ) > 0, A < ct'(s) < 
A Vs > 2 s 2, where A > A > 0 are constants. 

(hi) Let s^ G [0,si) and S 2 G (s 2 ,oo) denote the unique numbers with 
(t(s|) = (t(s2) and ( 7 ( 32 ) = cr(si), respectively. 


Note that the cases ct(s 2 ) > 0 and ( 7 ( 32 ) = 0 correspond to Figures]^ and 
respectively. In addition to Hypothesis (H), if we suppose cr(s) > 0 for all 
s G (si, S 2 ), then A(p) • p > 0 for all p G M”", and so (2.9) is satisfied by any 
Lipschitz solution rt to (1.2); but we do not explicitly assume such positivity 
for Case II unless otherwise stated. 

With Hypothesis (H), for each r G (( 1 ( 52 ),(tCsi)), let s_(r) G and 

s+(r) G ( 52 , 52 ) denote the unique numbers with r = (T(5±(r)). 
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The following theorem is the second main result of this paper that gen¬ 
eralizes those of iznE] to any dimension n > 1. 

Theorem 3.5. Let Mq = ||L)mo||l°°(o); -^o ~ ■si}, and\DuQ{xQ)\ G 

(s|,S 2 ) for some xq G O. Then for each r G (< 7 ( 52 ), ct(Mq)), there exists a 
number I = Ir £ (cr(s 2 ),r) sueh that for all r G {l,r) and all but countably 
many fi G {a{s 2 ),r), 73 G (r, (t(si)), there exist three disjoint open sets 
Oy, rij, C LIt with = \LIt\ and infinitely many Lipschitz 

solutions u to satisfying 

u G U Oy), ut = div{A{Du)) pointwise in U fly, 

\Du{x,t)\ < s-(n) y{x,t) £ ^T) (9n^, 

\Du{x,t)\ > s+(f 3 ) V(x,t) G rio® C 

|5| + |L| = |f7||, |5|>0, |L|>0, 


where 


and 


rigi = {(x,0) I X G n, \Duo{x)\ < s_(fi)}, 
iTff = {(x,0) I X G n, \Duo{x)\ > s+(f 3 )}, 

S = {(x,t) G I s-(ri) < |-Dtt(x,t)| < s_(r)}, 

L = {(x,t) G I s+(r) < |-Dn(x,t)| < s+(f 3 )}. 


In regard to this theorem, an explanation similar to Remark 3.2 can be 
made; but we omit this. As a byproduct, we also have the following simple 


existence result whose proof appears after that of Theorem 3.5 

Corollary 3.6. For any initial datum uq £ with Duq • n|gQ = 0, 

(1.2) has at least one Lipsehitz solution. 


We now restate Theorem 1.2 only for Case II. 

Theorem 3.7. Let 12 be convex and Mq = ||L)uo||L°°(n)- Assume further 
that f £ C'^([0, sf) U (^2 , 00 )). If Mq < si, then ( |l.2[ ) has a unique solution 
u £ satisfying \\Du\\loo(^q^) = ||T)no||Loo(n). 

Unlike Theorem |3 . 5| and its corollary, the convexity assumption on 12 can¬ 


not be dropped in this theorem. If < Mq < si, Theorems 3.5 and 3.7 


apply to give infinitely many Lipschitz solutions and one special Lipschitz 
(classical) solution to problem (1.2) under the same initial datum uq, re¬ 
spectively. However, a subtlety arises when Mq < Is there a Lipschitz 
solution to (1.2) other than the classical one in this case? Let us discuss 
more on this in the remark below. 

Remark 3.8 (Breakdown of uniqueness). In this remark, let a{s) satisfy 
(3.2) (t(s) > 0 VsG(si,S2] 

in addition to Hypothesi s (H ) (as in Figure]^, and let / be as in Theorem 


3.7 Then by Corollary |3.6[ for any uq £ C'^+"(12) with Duq ■ n|aQ = 0, 
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problem ( |1.2| ) has at least one Lipschitz solution. In particular, if Mq = 0 
with Mq := ||Z)uo||l°°(d)j that is, uq is constant in n, then by Proposition 
2.4, the constant function u = uq in Qt is a unique solution to (1.2). So a 


natural question is to ask if there is a number Mq > 0 such that (1.2) has a 
unique Lipschitz solution for all initial data uq G with DuQ-Ta\QQ, = 

0 and Mq < Mq. To address this question, define 


Mq = sup < Mq > 0 


Vuo G with Duq • n|gn = 0 and 

Mq < Mq, (1.2) has a unique Lipschitz solution 


which we may call as the critical threshold of uniqueness (CTU) for (1.2). 
The CTU Mq may depend not only on the profile a{s) but also on the 


dimension n, final time T, and domain Q of problem (1.2). By Theorem 3.5 
we have an estimate: 

0 < Mq* < si, 

which is valid for all dimensions n > 1, final times T > 0, and domains 
n satisfying (1.6). Our question now boils down to the point of asking if 
Mq* > 0. 

When n = 1, the answer is affirmative; in this case, the CTU Mq is 
positive and depends only on the profile cr(s). This fact can be derived by 
using an a priori estimate on the spatial derivatives of Lipschitz solutions. 
From 


when n = 1, any Lipschitz solution u to (1.2) satisfies 


(3.3) 


\'^x\\L°°{nT) — 
t{s) 




where c := inf 5 g(Q oo) > 0, by Hypothesis (H) and (3.2). Note 


O'('Sl) ^ o'(g2) ^ 


•si 


S2 


that is, > csi. As (t( 0) = 0 and (t(s) is strictly increasing on s G [0, 

there is a unique number Mq G (0, sJ) with u(M^) = csi. Suppose now that 
Then by Theorem 


u 


'0llL°°(D) 


< 


M'. 


3.7 


there exists a classical solution 


G to (1.2) with ||ttx||L^(o.r) = ||iiollL°°(D) < On the 


other hand, if u is any Lipschitz solution to (1.2), we have from (3.3) that 
< si- Modifying the profile cr(s) to the right of the threshold 
max{||uQ||£,cx3(Q), ||u 3 ;||L«)(Oy)} < si, both u and u become Lipschitz solutions 


to some uniformly parabolic problem of type (1.2) with monotone flux under 


the same initial datum uq. By Proposition 2.3, we thus have u = u in 
hence, the Lipschitz solution u to (1.2) is unique. Therefore, we have an 


improved estimate for the CTU Mq in n = 1: 


0 < M^ < Mq* < si- 


As a toy example for n = 1, we consider the general piecewise linear 
profile (t(s) of the Hollig type. Let si, S 2 , ki, k 2 and k^ be positive numbers 
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with S 2 > Si and —^ 2(52 — si) + fcisi > 0. Let us take 


(t(s) 


kis, 

-k 2 {s - si) + kisi, 

ksis - S 2 ) - k2{s2 - si) + kisi, 


0 < s < si, 
Si < s < S2, 
S > S 2 ; 


then, if —k2{s2 — si) + kiSi < ^332, we have 

0 < -si^ + ^(1 + ^) = M' < Mo* < s^ = 

kl S2 Kl 


k2 . k2. 


or if —k 2 {s 2 — si) + kiSi > /C 3 S 2 , it follows that 


0 < si^ = M^ < Mo* < s^ = -S2^ 



Unfortunately, estimate 
mensions n > 2. However, we still expect that the CTU Mq > 0 even in 
dimensions n > 2 for all T > 0 as long as H is convex. Since this study is 
beyond the scope of this paper, we leave it to the interested readers. 


(3.3) from [22] is not directly generalized to di- 


4. Some useful results 


This section prepares some essential ingredients for the proofs of existence 


theorems. Theorems |3.1| and 3.5 


4.1. Uniformly parabolic equations. We refer to the standard references 
(e.g., [53130]) for some notations concerning functions and domains of class 
(jk+a integer k > 0. 

Assume / G 00 )) is a function satisfying 

(4.1) 0 < f{s)+ 2sf'{s) <e Vs>0, 


where 0 > 0 > 0 are constants. This condition is equivalent to 0 < 
(s/(s^))' < 0 for all s G M; hence, 9 < f{s) < 0 for all s > 0. Let 

Mp) = /(|pP)p {p e 1^")- 

Then we have 

^lyip) = f{\p\^)^ij + 2/'(|p|^)piPj (i,j = 1,2, • • • ,n; p G M."-) 
and hence the uniform elUpticity condition'. 

n 

(4.2) e\q\‘^ < Al^{p)qiqj < Q\q\'^ ^ P, q ^ M”. 

*j=i 


Theorem 4.1. 

problem 


Assume (1.6). 


Then the initial-Neumann boundary value 


ut = div(A(Utt)) 
< du/dn = 0 
u{x, 0) = uo{x) 


in 

on dQ X (0,T), 
for X £ Q 


(4.3) 
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has a unique solution u G Moreover, if f G C^([0,(X))) and 

0 is convex, then the gradient maximum principle holds: 


(4.4) 


I I-Dull 1,00 — IIDUoIIloo(J^). 


Proof. 1. As problem (4.3) is uniformly parabolic by (4.2), the existence 
of unique classical solution u in follows from the standard 

theory; see |3Q1 Theorem 13.24]. 

2. To prove the gradient maximum principle (4.4), assume hereafter 

/ G C'^([0,oo)) and fl is convex. Note that, since A G a standard 

bootstrap argument based on the regularity theory of linear parabolic equa¬ 
tions |29[ l30] shows that the solution u has all continuous partial derivatives 
UxiXjXk and Uxit within Qt for 1 <i,j,k < n. 

3. Let V = \Du\‘^. Then, within Qt, we compute 

Av = 2Du ■ D{Au) + 2|T)^up, 
ut = div(A(Hu)) = div{f{v)Du) = f'{v)Dv ■ Du -|- f{v)Au, 

Dut = f"{v){Dv ■ Du)Dv + f' {v){D‘^u)Dv 
+f{v){D‘^v)Du + f'{v){Au)Dv -|- f{v)D{Au). 

Putting these equations into vt = 2Du ■ Dut, we obtain 

(4.5) Vt-C{v) - B ■ Dv =-2f{\Du\'^)\D'^u\‘^ <0 in Or, 
where operator C{v) and coefficient B are defined by 

C{v) = f{\Du\‘^)Av + 2f{\Du\‘^)Du ■ {D^v)Du, 

B = 2f{v){Dv • Du)Du + 2f{v){D^u)Du + 2f{v){Au)Du. 

We write C{v) = Y17j=i with coefficients atj = aij{x,t) given by 

aij = Ap^{Du) = f{\Du\‘^)6ij + 2f'{\Du\‘^)uxiUxj (f,j = l,-- - ,n). 

Note that on Or all eigenvalues of the matrix {atj) lie in [0, 0]. 

4. We show 

max v{x,t) = m.axv{x,0), 

which proves (4.4). We prove this by contradiction. Suppose 

(4.6) M := max v{x,t) > maxv{x,0). 

Let {xo,to) G Or be such that v{xo,to) = M; then to > 0. If xq G 0, 
then the strong maximum principle applied to (4.5) would imply that v is 
constant on O^p, which yields u(x,0) = M on O, a contradiction to (4.6). 
Consequently xq G 90 and thus u(xo,to) = M > v{x,t) for all {x,t) G Or- 
We can then apply Hopf’s Lemma for parabolic equations [36] to (4.5) to 


deduce dv{xo,tQ)/dn > 0. However, mainly due the convexity of 0, a result 
of [H Lemma 2.1] (see also [231 Theorem 2]) asserts that dv/dm < 0 on 
90 X [0,T], which gives a desired contradiction. □ 
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4.2. Modification of the profile functions. The following elementary 
results can be proved in a similar way as in [HI Ha]; we omit the proofs. 

Lemma 4.2 (Case I: Perona-Malik type; see Figure]^. Assume Hypothesis 
(PM). For every 0 < n < r 2 < (t{so), there exists a function a G C'^([0, oo)) 
such that 


(4.7) 


a{s) = cr(s), 

< a{s) < a{s), 

6 < a'{s) < 0, 


0 < s < s_(ri), 
s_(ri) < s < s+(r 2 ), 
0 < s < oo 


for some constants @ > 6 > 0. With such a function a, define f{s) = 
a{^)/^/s (s > 0) and /(O) = /(O); then f G C'^([0,oo)) fulfills condition 

dtll). 



Figure 6. Case I: Perona-Malik type of prohle a{s) and 
modified function 5'(s). 


Lemma 4.3 (Case II: Hollig type; see Figure]^. Assume Hypothesis (H). 
For every (t{s 2 ) < ri < r 2 < cr{si), there exists a function a G C^+“([0, oo)) 
such that 


(4.8) 


/ 

d(s) = fT(s), 
a-(s) < (t(s), 

< a(s) > (t(s), 
ci(s) = cr(s), 

9 < a'(s) < 0, 


0 < s < s_(ri), 
s_(ri) < s < s_(r 2 ), 
s+(ri) < s < s+(r 2 ), 
s+{r 2 ) < s < oo, 

0 < s < oo 


for some constants Q > 6 > 0. With such a function a, define f{s) = 
^(\/ s)/\/s (s > 0) and /(O) = /(O); then f G C'^'''"([0, oo)) fulfills condition 
(4.1). Moreover, if f € (^^([Ojsf) U {s^ , oo)) in addition, then a, f can be 
chosen to he also in C'^([0,oo)). 
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Figure 7. Case II: Hollig type of profile (t(s) and modified 
function ct(s). 


4.3. Proofs of Theorems 3.4 and 3.7 


Let Mq = ||L)no||L°°(r 2 )- In Case 
I (Perona-Malik type), we nave 0 < Mq < sq; so we can select 0 < ri < 
T 2 < t('So) such that S-(ri) = Mq. In Case II (Hollig type), we have 0 < 
Mq < si; so we can select ct(s 2 ) < n < r 2 < (t(si) such that s_(ri) > Mq. 


Use such a choice of ri, r 2 in Lemma 4.2 (Case I), Lemma 4.3 (Case II 


to obtain a functio n / as stated in the lemma. Let A{p) = f{\p \'^)p. 


For t his A(p), problem (4.3) is uniformly parabolic; hence, by Theorem 
(4.3) has a unique solution u G satisfying 


4.1 


Since A{p) = A{p) on |p| < s_(ri) and S-(ri) > Mq, it follows that 
A{Du{x,t)) = A{Du{x,t)) in Qt] this proves rt is a classical solution to 
problem (1.2). On the other hand, we easily see that any classical solution 


to (1.2) satisfying \\Du\\ioo^q^) = Mq is also a classical solution to (4.3) and 
hence must be unique. This completes the proofs of Theorems 3.4 and 3.7 


4.4. Right inverse of the divergence operator. To deal with linear 
constraint divn = u, we follow an argument of O Lemma 4] to construct 
a right inverse TZ of the divergence operator: div77 = Id (in the sense 
of distributions in Ht). For the purpose of this paper, the construction 
of IZ is restricted to box domains, by which we mean domains given by 
Q = Ji X J2 X ■ ■ ■ X Jn, where Jj = (o^, 6*) C M is a finite open interval. 

Given a box Q, we define a linear operator IZn- L°°{Q) —)■ L°°{Q]W^) 
inductively on dimension n as follows. If n = 1, for u G we define 

V = IZiu by 

nxi 

v{xi) = / u{s)ds {xi G Ji). 

J ai 
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Assume n = 2. Let u G x J 2 ). Set u{xi) = u{xi, s) ds for xi G Ji. 

Then u G L°°(Ji). Let v = that is, 

r-xi rxi j-b2 

h(a:i) = / u{s)ds = / / u{s,T)dTds (xi G Ji). 

^ai J ai J a2 

Let p 2 G C^{a 2 ,b 2 ) be such that 0 < /92(s) < /® 2 ('S)(is = 

1. Define v = TZ 2 U G L°°{Ji x J 2 ;M^) by u = with u^(xi,X 2 ) = 

P 2 {x 2 )v{xi) and 

rX2 rX2 

V^{xi,X2)= / tt(xi, s)(is — -u(xi) / P 2 {s)ds. 

J a2 J a2 

Note that if u G x J 2 ) then u G hence u = TZ 2 U G 

X J 2 ; M^) and div v = u a.e. in Ji x J 2 . Moreover, if n G C'^(Ji x J 2 ) 
then V is in C^{Ji x J 2 ;M^). 

Assume that we have defined the operator TZn-i- Let u G L°°{Q) with 
Q = Ji X J 2 X ■■■ X Jn and x = {x', Xn) G Q, where x' ^ Q' = Jix ■ ■ ■ x Jn-i 
and Xn G Jn- Set u{x') = u{x', s) ds for x' G Q'. Then u G L°°{Q'). 
By the assumption, v = TZn-iu G L°°{Q 'is defined. Write v{x') = 
{Z^{x'), ■ ■ ■ , Z^~^{x')), and let pn G C^{an,bn) be a function satisfying 
0 < Pnis) < j, and pn{s)ds = 1. Define v = TZnU G L°^{Q‘,W^) as 
follows. For X = (x', x^) G Q, u(x) = (u^(x), u^(x), • • • , u”'(x)) is dehned by 

v^{x', Xn) = Pn{xn)Z^{x') (A: = 1, 2, • • • , n - 1), 

rxn rxn 

v"'{x',Xn)= / u{x', s)ds — u{x') / Pn{s)ds. 

J Clji J Cl-f^ 

Then TZn- L°^{Q) —)■ L'^ {Q]W^) is a well-defined linear operator; moreover, 

(4.9) ||7^„u||^oo(g) < Cn (|«/l| I Jn|)||u||ioo(Q), 

where (7,1 > 0 is a constant depending only on n. 

As in the case n = 2, we see that if u G W^’^{Q) then v = TZnU G 
M”) and divu = u a.e. in Q. Also, if n G C^{Q) then v = TZnU is 
in (7^((5;M"'). Moreover, if tt G W^’°°{Q) satisfies fQu(x)dx = 0, then one 

can easily show that v = TZnU G PFq’°°((5; 

Let / be a finite open interval in M. We now extend the operator TZn to 
an operator TZ on L°°{Q x I) by defining, for a.e. (x,t) £ Q x I, 

(4.10) (7^u)(x,^) = (7^„u(•,^))(x) y u£ L°^{Q X I). 

Then TZ: L°^{Q x /) —)■ L°^{Q x I;M"') is a bounded linear operator. 

We have the following result. 

Theorem 4.4. Let u G Wq’°°{Q x I) satisfy fg u(x, t) dx = 0 for all t £ I. 
Then v = TZu £ x I;W^), divu = u a.e. in Q x I, and 

(4-11) II'^i||l°°(Qx/) < Cn (|7l| -I- • • • -I- \Jn\)\\ut\\L°°{QxI), 



FORWARD-BACKWARD PARABOLIC EQUATIONS 


23 


where Q = Ji x ■ ■ ■ x Jn and Cn is the same constant as in (4.9). Moreover, 
ifuG C\Q^) then v = nu£ 

Proof. Given u G Wq’'^{Q x /), let p = TZu. We easily verify that v is 
Lipschitz continuous in t and hence vt exists. It also follows that vt = TZ{ut). 
Clearly, if u(x, t)dx = 0 then v{x, t) = 0 whenever t G dl or x £ dQ. 

This proves v G x and the estimate (4.11) follows from 

(4.9). Finally, from the definition of TZu, we see that if G C^(Q x I) then 

V = TZuGC\Q^;R^). □ 

5. Geometry of the relevant matrix sets 


Let A{p) be given by (1.3). We assume Hypothesis (PM) or (H) unless 


one is chosen specifically. Recall the definition (|2.3|) with 1 = 0: 

trR = 0 


(5.1) Ko = K{0) = 


p c 
B A{p) 


p G M”, c G M, R G 


Under Hypothesis (PM) or (H), certain structures of the set Kq turn 
out to be still quite useful, especially when it comes to the relaxation of 
homogeneous partial differential inclusion Vuj{z) G Kq with z = (x, t) and 
a; = ((p,V’)- We investigate these structures and establish such relaxation 
results throughout this section. 

5.1. Geometry of the matrix set Kq. We study some subsets of Kq, 
depending on the different types of profiles. 

Case I; Hypothesis (PM). (See Figure]^) In this case, we assume the 
following. Fix any two numbers 0 < ri < r 2 < ( 7 ( 50 ), and let Fq = ^^^^^ 2 ( 0 ) 
be the subset of Kq defined by 

p c \ p G IpI G (s_(ri),s_(r 2 )) U (s+(r 2 ),s+(ri)), 

B A{p)j c G M, R G trR = 0 

We decompose the set Fq into two disjoint subsets as follows: 

p c \ p G M*", IpI G (s_(ri),s_(r 2 )), 

R A{p)) c G M, R G tr R = 0 


Rn = 


R_ = 


F. = 


p c 
B A{p) 


p G W^, IpI G {s+{r 2 ),s+{ri)), 
c G M, R G trR = 0 


Case II: Hypothesis (H). (See Figure]^) In this case, we assume the 
following. Fix any two numbers ct(s 2 ) < ri < r 2 < cr(si), and let Fq = 
Fr 


ri,r2(0) be the subset of Kq given by 


Rn = 


p c 

R H(p) 


R_ = 


p G M”, IpI G (s_(ri),s_(r 2 )) U (s+(ri), s+(r 2 )), 
c G M, R G tr R = 0 

The set Rq is also decomposed into two disjoint subsets as follows: 

p c \ p G IpI G (s_(ri),s_(r 2 )), 

R A{p)J c G M, R G trR = 0 
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p G W^, \p\ G (s+(ri),s+(r 2 )), 1 
c G M, -B G trB = 0 J ' 

In order to study the homogeneous differential inclusion Vuj{z) G Kq, we 
hrst scrutinize the rank-one structure of the set Fq C Kq. We introduce the 
following notation. 

Definition 5.1. For a given set E C L{E) is defined to be the 

set of all matrices ^ G that are not in E but are representable 

by ^ = A^i -|- (1 — A )|2 for some A G (0,1) and ^2 G B with rank(^i —^ 2 ) = 
1 , or equivalently, 

L{E) = ^ E \ ^ + t±ri G E for some t_ < 0 < and rank?] = 1}. 

For the matrix set Fq, we define 

R{Fo)= U 

C± eE±, rank(^+ )=1 

where is the open line segment in joining Then 

from careful analyses, one can actually deduce 

(5.2) L(Bo) = R{Fo) U L{F+) in Case I 

and 



(5.3) 


L(Bo) = B(Bo) in Case II. 


In (5.2), due to the backward nature of cj(s) on (s+(r 2 ), s+(ri)) for Case 
I, the set B(F+) turns out to be non-empty. On the other hand, as only 
forward parts of a are involved in Fq for Case II, no such set appears in 
(5.3). Regardless of this discrepancy, we can only stick to the analysis of the 
set R{Fq) for both cases towards the existence results. Theorems 3.1 and 


We perform the step-by-step analysis of the set R{Fq) for both cases 
simultaneously. 


5.1.1. Alternate expression for R{Fq). We investigate more specific cri¬ 
teria for matrices in R{Fo). 

Lemma 5.2. Let ^ G Then ^ G R{Fq) if and only if there 

exist numbers t- < 0 < t+ and vectors q, ^ £ ME with |(]| = 1 , 7 • g = 0 such 

that for each 5 G M \ {0}, ^ ~ ^ ^ ^ 


p 

B 


Proof. Assume | = 
t_ < 0 < and 7 is a rank-one matrix given by 


G R{Fo). By definition, ^ -|- t±fi G F±, where 


7 ] = 


a 


® (g,^) = 


aq ab 
a® q ba 


a^ F lap / 0, b^ + |gp / 0, 
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for some a, 6 G M and a, q ^ M”; here a® q denotes the rank-one or zero 
matrix {aiqj) in Condition ^ -|- t±fj G F± with < 0 < is 

equivalent to the following: For Case I, 


(5.4) 


tr i? = 0, a ■ q = 0, A{p + t±aq) = /3 + t±ba, 

\p + t+aq\ G (s+(r 2 ),s+(ri)), \p + t-aq\ G (s_(ri), s_(r 2 )). 
For Case II, 

tr B = 0, a ■ q = 0, A{p + t±aq) = {3 + t±ba, 


(5.5) 


\p + t+aq\ G (s+(ri),s+(r 2 )), \p + t-aq\ G (s_(ri), s_(r 2 )). 


Therefore, aq ^ 0. Upon rescaling fj and t±, we can assume a = 1 and 
|g| = 1 ; namely, 


T] = 


q b 
a® q ba 


\q\ = l, a-q = 0. 


We now let 7 = ba. Let 6 G M \ {0} and 

q b 


^ 'i7<»9 7, 


From (5.4) or (5.5), it follows that ^ + t±r/ G F±. 

The converse directly follows from the definition of R{Fq). 


□ 


5.1.2. Diagonal components of matrices in R{Fo). The following gives 
a description for the diagonal components of matrices in R{Fq). 


Lemma 5.3. 

(5.6) R{Fo) = 


p c 
B p 


c G M, 5 G 


trS = 0, (p,/3) G S 


for some set S = Sri,r 2 C 

Proof. Let {c,B), {c',B') G M x 
define 


be such that tr B = tr B' = 0, and 


SicB) = |(p,/3) e I G i?(Fo)| , 


S{c',B') — \ iP,P) e 


B' /3 


G RiFo) } . 


It is sufficient to show that S(c,B) = 'S{c',b') =■ <S. Let (p,/3) G S(c,b)j that 
^ G R{Fo). Then C± •= ? + t±V G F± for some < 0 < 

—t — 


is, ^ = 

and ra: 
and that 


and rankr/ = 1. Observe that f = A^+ -|- (1 — A)^_ with A = -j—p- G (0,1) 


e' := 


p 


B' 13 


= C + 


0 c 
B 0 


= Ae+ + (1 - A)|_ 
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where c = d — c, B = B' — B, and 


C± + • Since £ 


F± and tiB = 0, we have G F±, and so G R{Fq). This implies 
{p,/3) G hence 5(c,b) C S(^c',b')- Likewise, S(^c',B') F S^^b)', that is, 

S{c,B) = S{^c',B')- n 


5.1.3. Selection of approximate collinear rank-one connections for 

R{Fq). We begin with a 2-dimensional description for the rank-one connec¬ 
tions of diagonal components of matrices in R{Fq) in a general form. The 
following lemma deals with Cases I and II in a parallel manner. 


Lemma 5 . 4 . For all positive numbers a, b, c with b > a, there exists a 
continuous function 

hi{a, b, c; •, •, •) : = [0, a) x [0, oo) x [0, c) —> [0, oo) (Case I) 

/i 2 (a, b, c; •, •, •) : Ia^b,c = [Oj ^ [0, 6 — a) x [0, c) —> [0, oo) (Case II) 

with hi{a, b, c; 0 , 0 , 0 ) = 0 satisfying the following: 

Let (5i, 82 and rj be any positive numbers with 

0 <a — Si<a<b<b + 62, 0 < c — rj < c, (Case I) 

0<a — (5i<a<6 — ( 52 < 6 , 0 < c — p < c, (Case II) 

and let Ri G [a —( 5 i,a], R2 G [8,8 + 62] (Case I), R2 G [8 — 62,6] (Case II), 
and Ri, R2 £ [c — p,c]. Suppose 6 £ [—7r/2,7r/2] and 

(-Ri(cos(| -h 6 '),sin(| -ht*)) -.R 2 (cos(| - 6 '),sin(| - 6 )))) 

• (^i?i(cos(| -h 6 '),sin(| -ht*)) - ii 2 (cos(| - 6 '),sin(| - 6 *))^ = 0 . 
Then — R2, and 

max||(0,a)-i?i(cos(|-k(9),sin(^-F6»))|, | ( 0 , 6 )-i ?2 (cos(|- 6 '), sin(|- 6 l)) | 

I ( 0 , c) - .Ri (cos(^ -h 6 »), sin(| + 6 *)) |, | ( 0 , c) - .R 2 (cos(^ - 6 »), sin(| - 6 »)) 11 

< hi{a, b, c; 61,62, p), (Case I) 

max||( 0 ,a)-i?i(cos(|-k 6 '),sin(|-F 6 »))|, | ( 0 , 6 )-i ?2 (cos(|-6'), sin(|-6l)) | 

I (0, c) - .Ri ( cos(^ -h 6»), sin(| + 6))) |, | (0, c) - .R 2 ( cos(^ - 6»), sin(| - 6))) 11 

< /i 2 (a, b, c; 6 i, 62 ,p). (Case II) 

Proof. By assumption, 

0 = {Ri{— sin 0 , cos 6 ) — R2 (sin 9 , cos 9 )) ■ {Ri{— sin 0 , cos 9 ) — R2 (sin 9 , cos 9 )) 
= {—{Ri + R2) sin0, (i?i — R2) cos 9 ) ■ {—{Ri + R2) sin0, (Ri — R2) cos 9 ) 
= [Ri + R2 ){Ri + R2) sin^ 9 + {Ri - R2){Ri - R2) cos^ 9 , 
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that is, 


(i ?2 - Ri){Ri - R 2 ) cos^ 6 = {Ri + R 2 ){Ri + R 2 ) sin^ 9; 


hence, 6 / ±|, Ri > R 2 , and 


0 = ± tan 


So 


\9\ < tan 


/ {b-a + 5i + 52 )rj 


( 1 {R2-Ri){Ri-R2) \ 

(i?l + R 2 ){Ri + R 2 ) J 

=: gi(a,b,c; 61 , 62 ,-rj), (Case I) 


\ 6 \ < tan 


-1 


2 {a + b- (5i)(c - rf) 
(6 - a + 61 ) 7 ] 


=: g 2 {a,b,c; 6 i, 62 ,r]). (Case II) 


^ y 2 {a + b - 61 - 62 ){c - g) ^ 

Note that the function gi{a,b,c; ■,■,■) '■ lab [0)^/2) is well-defined and 
continuous and that gi{a,b,c-, 61,62, g) = 0 for all {6i,62,g) & la be with 


g = 0 . 

Observe now that 


1(0, a) - i?i(cos(| (9),sin(| -h (9))| 


< max{|( 0 , a) — a(— sin 9 , cos 0 )|, | ( 0 , a) — (a — 6i)(— sin 9 , cos 6 *)|} 

= max I sin^ 9 + a^(l — cos 9)^, \J {a — 61)“^ sin^ 9 + {a — {a — (5i) cos0)2| 

= max |\/2a\/l — cos 9 , 

< max|\/2av^l - cos{gi{a,b,c;6i,62,g)), 

\/{a- (5i)2 a 2 - 2a{a - (5i) cos{gi{a,b, c; 6i,62,g))'^ 

=: hi^i{a,b,c; 61,62, g), (i=l, 2) 

1 ( 0 , 6 ) -ii 2 (cos(^ - 0 ),sin(^ - 0 ))| 

< max{ 1(0,6) — 6 (sin 9 , cos 0 ) | , | (0, 6 ) — {b + 62 ) (sin 9 , cos 9 ) \} 

= max I y^ 62 sin^ 0 -|- 6^(1 — cos0)^, y/(6 -|- 62)^ sin^ 6* -|- (6 — (6 -|- 62) cos 0)^1 
= max I \f 2 b\/ 1 — cos 9 , \/(b^ 62Y ^b^ — 26(6 -|- 62) cos 6*| 
max |\/26v^l - cos{gi{a,b, c; 6i,62,g)), 

a /(6 -I- 62)2 -I- 62 - 26(6 -I- 62) 003(31(0, 6 , c; 61,62, 

=: 6-1, 2 ( 0 , 6, c; 61, 62, g), (Case I) 

1(0,6) -ii 2 (cos(^ -0),sin(^ -0))| 


< 
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< max 


|\/26v^l - cos{g2{a,b, c; 6i,52,r])), 

\/{b- 62)'^ + 62 - 26(6 - 62)cos{g2{a,b,c- (5i, <52, i/))| 

=: /i 2 , 2 (a, b, c; 61,62,1]), (Case II) 


|(0,c) -«i(cos(^ + (9),sin(| + (9))| 


< 


max |\/ 2 c-v/l - cos{gi{a,b,c; 61,62,1])), 


Vic- viY + c 2 - 2 c(c - ry) cos{gi{a,b,c-, (5i, < 52 , r/))| 

=: hi^2.ia,b,c-, 61,62,V, ( 1 = 152 ) 

1 ( 0 ,c) -^ 2 (cos(^ - 6 '),sin(| - 6 »))| < hi^3{a,b,c-,6i,62,v)- ( 1 = 1 , 2 ) 

Define hi{a,b,c-, 61,62, tj) = maxi<j< 3 /ijj(a, 5, c; (5i, ^ 2 ) 1 ?) for i = 1,2 cor¬ 
responding to Cases I and II, respectively. Then it is easy to see that 
the function hi{a, b, c; •,•,•) ■ Vhc ^ i® well-defined and satisfies the 

desired properties. □ 

We now apply the previous lemma to choose approximate collinear rank- 
one connections for the diagonal components of matrices in R{Fq). 


Theorem 5 . 5 . Let p± G M” satisfy 

s-{ri) < \p_\ < s-(r2) < s+(r2) < |p+| < s+(ri), (Case I) 
s_(ri) < |p_| < s_(r2) < s+(ri) < \p+\ < s+(r2), (Case II) 
and {A{p^) — A{p-)) ■ (p+ — p-) = 0. Then there exists a vector V G 
such that, with p^ = s±(r2)C°; ^(p±) = 'C2V, we have 

max{|p° - p_\,\pI - p+\,\A{pV - A{p_)\,\A{pV - A{p+)\} 

< hi{s-{r2),s+{r2),r2-,s-{r2) - s_ (n), s+(ri) -s+(r2),r2 -n), (Case I) 

max{|p° - p_\,\pI - p+\,\A{pV - A{p_)\,\A{pV - A{p+)\} 

< /I2(s-(r2),s+(r2),r2;s_(r2)-s_(ri),s+(r2)-s+(ri),r2-ri), (Case II) 
where hi, /12 are the functions in Lemma \ 5 . 4 \ 


Proof. Let S 2 denote the 2-dimensional linear subspace of M” spanned by 
the two vectors p±. (In the case that p± are collinear, we choose S 2 to be 
any 2-dimensional space containing p±.) Set 




p+ I p- 

b+l b-l 

I p+ I p- I 

I b+l b-M 


GS"“^nS2. 


Since the vectors p±, A{p±) and all lie in S2, we can recast the problem 
into the setting of the previous lemma via one of the two linear isomorphisms 
of S 2 onto with correspondence -H- (0,1) G M^. Then the results 
follow with the following choices in applying Lemma |5.4[ a = S-{r 2 ), b = 
s+{r2), c = r2, ( 5 i = s_(r2) - s_(ri), 62 = s+(n) - s+(?’2) (Case I), 62 = 
s+{r2) - s+{ri) (Case II), r] = r2 -ri, Ri = \p_\, R2 = \p+\, Ri = o'(IP-I), 
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-^2 = crdP+Di and 6 G [0,7r/2] is the half of the angle between and 

p_. □ 


5.1.4. Final characterization of R{Fq). We are now ready to establish 
the result concerning the essential structures of R{Fq). For this purpose, it 
suffices to stick only to the diagonal components of matrices in R(Fo). 

Theorem 5 . 6 . Let 0 < r 2 < cr(so) (Case I), it(s2) < r 2 < (Case 

II). Then there exists a number I 2 = lr 2 £ (0,r2) (Case I), a number 
I 2 = lr 2 £ (< 7 (^ 2 ), F 2 ) (Case II) such that for any I 2 < ri < r 2 , the set 
S = S, 


ri,r2 


c in (|5.6l) satisfies the following: 


(i) sup(p_^)g 5 IpI < s+(ri) (Case I), sup(p_^)g 5 |p| < s+(r 2 ) (Case II) 
and sup(p^^)g 5 |/3| < r 2 ; hence S is bounded. 

(ii) S is open. 

(iii) For each {po,f3o) G S, there exist an open set V CC S containing 


(po,/?o) o.nd functions g : V —)• S"" 7 : V 


t± : V 


with 7 
p c 
B fi 


g = 0 and t- < 0 < t+ on V such that for every ^ = 
G R{Fq) = R{Fri^r2i0)) with {p, f3) G V, we have 


^ + t±p € F±, 


where t± = t±{p,(3), p = 
is arbitrary. 


q{p,fi) b 


and 5/0 


Proof. Fix any 0 < r 2 < (j{so) (Case I), 17 ( 52 ) <r 2 < cr(si) (Case II). For 
the moment, we let ri be any number in (0, r 2 ) (Case I), in ( 17 ( 52 ), ^ 2 ) (Case 
II) and prove (i). Then we choose later a lower bound I 2 = lr 2 £ ( 0 , 72 ) 
(Case I), I 2 = lr 2 £ (< 7 ( 52 ), r 2 ) (Case II) of ri for the validity of (ii) and 
(iii) above. 

We divide the proof into several steps. 

1. To show that (i) holds, choose any (p,/3) G S. By Lemma 5.3, ^ := 

p 0 

O /3 

'p± 

B± ^{\P±\)^\ 

and a number 0 < A < 1 such that ^ = A^+ + (1 ~ So 


G R{Fq) = i?(Fr^_r2(0))) where O is the nx n zero matrix. By the 

c± 


definition of R{Fo), there exist two matrices = 




\p\ = \>^P+ + (1 - A)p-| < s+(ri), (Case I) 


IpI = |Ap+ + (1 - A)p_| < 5 +( 72 ), (Case II) 


l/3| 


At^(|p+I)j^ + (1 - A)(7 (|p_|) 


< 72 ; 


hence, sup(p^^)g 5 |p| < 5+(ri) (Case I), sup(p_^)g 5 |p| < 5+(r2) (Case II), 
sup(p^/ 3 )g 5 |/3| < 72, and S is bounded. So (i) is proved. 
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2. We now turn to the remaining assertions that the set S = Sri,r 2 fulfills 
(ii) and (iii) for all ri < r 2 sufficiently close to r 2 . In this step, we still 
assume ri is any fixed number in (0,r 2 ) (Case I), in {cr{s 2 ),r 2 ) (Case II). 

Po 0 

O P 02 

5.2| that there exist numbers sq < 0 < to and vectors go, Jo £ with 


Let {po,Po) G S. Since ^0 := 


G R{Fo), it follows from Lemma 


Igol = 1, 7 o • = 0 such that .^0 + sqtjo G F- and ^0 + toVo £ where 

1 ® and 6 7^ 0 is any fixed number. Let q'r, = tnon 7^ 0, 

70 Joj ^ HO OHO ^ , 

7o = ^o7o, and Sq = so/to < 0; then 


Vo = 


(5.7) 


To • ^0 = s-{ri) < \po + Sq^oI < s_(r 2 ), 
s+(r2) < \po + q'o\ < s+in), (Case I) 
s+(ri) < \po + q'qI < s+(r2), (Case II) 

^i\Po + gQgoD ipo+^l^gil = /^o + ^o7o, 

^^{\Po + q'o\)§^\=^o + io- 

Observe also 

fo - So > |(Po + io?o)| - \{Po + sogo)| > s+(r2) - s_(r2), 
to - So > |(po + io?o)| - |(Po + sogo)| > s+(ri) - s_(r2). 
Next, consider the function F defined by 

/a(|p + s'g'|)|±f^ - 13-s'i 
<yi\p + <i'\)^\-P-i 


(5.8) 


(5.9) F{^',q',s';p,/ 3 ) = 


(Case I) 
(Case II) 


G 


V 


y-g' 


for all 7', q', p, (3 G M"' and s' G M with s_(ri) < \p + s'g'| < s_(r 2 ), 
s+(r 2 ) < \p + q'\ < s+(ri) (Case I), s+(ri) < \p + q'\ < s+(r 2 ) (Case 
II). Then F is in the described open subset of and the 

observation (5.7) yields that 

F{jo:%^So-,po,/3o) = 0. 


Suppose for the moment that the Jacobian matrix D^y^gi^gi'^F is invertible 
at the point {Jq, q'o, Sq]Po, Po)', then the Implicit Function Theorem implies 
the following: There exist a bounded domain V = C containing 

{Po, Po) and functions g, 7 G M"', s G M of (p, /?) G V such that 

l{po,Po) = 7 o> g(Po, Po) = q'o, s(po, Po) = s'o 


s(p,P) < 0, s_(ri) < jpF s(p,p)q(p,P)l < s_(r 2 ), 
s+(r 2 ) < jp + q(p,P)j < s+(ri), (Case I) 
s+(ri) < \p + q{p,P)\ < s+(r 2 ), (Case II) 
F{jip, P),q{p, P),s{p, P);p, P) = 0 V(p,/3) G V. 


and that 
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Define functions 


7 = -pu, = t+ = \q\ in V; 

\q\ \q\ 

then 

s_(ri) < \p + t-q\ < s-(r 2 ), 
s+(r 2 ) < \p + t+q\ < s+(ri), (Case I) 
s+(ri) <\p + t+q\ < s+(r2), (Case II) 

a{\p + t±q\) ^ ^ =13 + t± 7 , |g| = 1, 7-9 = 0, t_ < 0 < t+, 

\p + t±q\ 

where {p,l3) G V, 7 = 7 (P,^), 9 = 9(p,/3), and = t±{p,l3). 

Let (p,/3) G V, 5 E tiB = 0, b,c £ R, b 0, q = q{p,P), 


7 = liP,/3), t± = t±{p,l3), C = 


p c 
B p 


, and 7 = ( 1 


b7<»9 7 


Then 


:= C + t± V G F±. By the definition of R{Fq), ^ G (^_,,^+) C R{Fq). By 
we thus have (p,/3) G 5; hence V C 5. This proves that S is 


5.3 


Lemma 

open. Choosing any open set V CC V with (po,/3o) G V, the assertion (iii) 
will hold. 

3. In this step, we continue Step 2 to deduce an equivalent condition for 
the invertibility of the Jacobian matrix Dj-y g/,,/)T at ( 79 , 9 o,'Sg;po,/^o)- By 
direct computation, 




D 




oT = 


Ml 


|y[j-(n-|-n-|-l) X (n-|-n-|-l) 


9 7 

where R is the n x n identity matrix. 


Mg, = s' |^o-'(|p + s'g'l) - 
= (^'i\P +s'q'l) - 


cT(|p + s'g'|)^ p + s'g' ^ p + s'g' ^ j a{\p + s'q'\) ^ 


Ip + s'q' 
o-(|p + s'g'l) 


|p + s'g'| lp + s'9'l 


Ip + s'g'l 


Ip + s'g'l y V Ip + ^'q' 


p + s'q' P +'S' 9 ^ _|_<7(|p + s'g'l) 


lp + s' 9 'l |p + s'g'| 


Here the prime only in cr' denotes the derivative. For notational simplicity, 
we write ( 7 ', g', s';p,/3) = ( 79 , gg,'S' q;po,/ 3o)- Applying suitable elementary 
row operations, as s' < 0 , we have 


F(-y',q',s') F 



Ms' 

Ml - FMs' 


^s' _ 


+ 'F{Ms'Y + --- + 'F{Ms'Y 


'-s'in Ms' 

O s'Ml - Ms' 

0 s' 7 ' + q'i{Ms'Y + • • • + q'n{Ms') 


^s'_ 

" q'-oj-', 


9 '± 7 '. 
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where O is the n x n zero matrix, and is the ith row of Mg/. Since 

\q'\ = to, 'y' ■ q' = 0, and s-(ri) < \p + s'q'\ < s_(r2), we have 

(^{\P +s'q'\) ^2 

\p + s'q'\ ° 


/ - f fn : / /\\ o-(|p + s'g'|)^^ f p + s'q' 


= tg I cos^ 0 O' {\p + s g I) + (1 — cos"^ 6 ) 


7| • g I + ' . :/' t, 


> 0 , 


\p + s'q'\ 

where 9 ' G [ 0 , vr] is the angle between p + s'q' and q'. Observe here that 
the forward part of a in the definition of F- becomes essential to guarantee 
that o'{\p + s'q'\) > 0 . After some elementary column operations to the last 
matrix from the above row operations, we obtain 

'-s'In Mg! - Ng! U]~i 
O s' AI^ — Alg/ -\- Ngf —CO / 


{-y',q',s') ^ 


0 0 


where the jth column of Ng/ G is 


s'y'.+q'-{M^,)j _ 




uj^,. So D(^Y^g/^gi-^F is 


invertible if and only if the n x n matrix Mi — ^Mgi + ^Ngi is invertible. 

We compute 

AY ^ AY I ^ hi I 'h + P + g' ^ P + l' 

o-dP + g'l). f ,n . / /|i o■(|p + s'g'l)^ p + s'g' ^ p + s'g' cj(|p + s'g 

j^In-[o {\p +s q\)- 


\p + q'\ 

1 


\p + s'q'\ J \p + s'q'\ |p + s'g'| Ip + s'g 






i + (o-'(|p + s'g I)- 


o'dp + 's'g'l)^ ( p + s'q 


= (ai - agi)In + {h - ai) 


Ip + s'q'l 
p + q' ^ p + q' 


Ip + s'q' 




p + s'q' 

Ip + ■s'g'l 


+ 


|p + g'| |p + g'| 


(65/ Ogf^ 


p + s q p + s q 1 -„ + 

]——77 ( 8 ) ]-—777 + -- - U ! , (g) UjJ ,, 

p + s'g' P + 'Sg g -w, 


and set (with an assumption ai / a^/) 

B :=-^-(Ml - ^Mg, + 

ai — Og' s s' 


—In + 


61 - ai p + q' p + q' 


ai-Os'lp + g'l |p + g'| 

bg! — Ogi p + s'q' 


p + s'q' 


+ 


where 


ai - Og' \p + s'q'\ Ip + s'g'd {ai - agi)q'■ u: 

o{\p + s'q'\) 


—oj , (g) ojZ, 


a.i = 


bg, = o'{\p + s'ci\)', 


Ip + s'q' 

then D^y^gi^gi-jF is invertible if and only if the matrix B G is invertible. 

4 . To close the arguments in Step 2 and thus to hnish the proof, we 
choose a suitable I2 = lr2 ^ (0)^2) (Case I), I2 = lr2 ^ (o(s2),r2) (Case 
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II), depending on r2, in such a way that for any ri G {l2,r2), the matrix 
B, determined through Steps 2 and 3 for any given {po,I 3 q) £ S = Sri,r2, is 
invertible. 

First, by Hypothesis (PM) or (H), r2 £ ( 0 ,r 2 ) (Case I), r2 £ (o'(s 2),F2) 
(Case II) can be chosen close enough to r2 so that 

a{k) a{l) 


< 


VZ G [s_(f2),S-(F2)], 


k I 

\/k £ [s+(r2),s+(f2)] (Case I), Vfe G [s+(f2), s+(r2)] (Case II). 
Then define a real-valued continuous function (to express the determinant 
of the matrix B from Step 3 ) 

, - -'(H)-tP 

DET(», t.,,, 7) = det ^ 


l“l 


hi 


u 

lttl 


o-jhl) _ o-(hl) |u| hi 


hi 


hi 








hi 


hi 


hi 


\ n tn 17 crd'^Dx/ V ^ V irdnl) 

^ ® ^ M •« M + ^ 

|t?| pI |z;| pI 

on the compact set A 4 of points {u, v, q, 7) G M" x M”" x x M” with 

IdI G [s_(f2),s_(r2)], ItI < 1 , 

|u| G [s+(r2),s+(f2)] (Case I), |u| G [s+(f2), s+(r2)] (Case II). 
With k = s+{r2) and I = s_(r2), for each q £ 

a'{k) - ^ ^ 

(T(fc) cr{l) 

“fc B 

since cr'{k) 7^ 0 and hence the fraction in front of q® q is different from —1. 
So 

d := min |DET(A;(7, Ig, g, 0)| > 0. 

Next, choose a number 5 > 0 such that for all {u, v, q, 7), {u, v, q,^) £ M. 
with |tt — m|, |ti — til, \q — q\, I7 — 7I < 6 , we have 


DET(fcg', Iq, q, 0 ) = det (In + 


9 7 ^ 0 , 


( 5 . 10 ) |DET(ri, n, q, 7) — DET(tt, v, ^,7)] < d/2. 

Let I2 £ {r2,r2) be sufficiently close to r2 so that for all ri G {12,^2), 
/ii(s_(r2),s+(r2),r2;s-(r2)-s-(ri),s+(ri)-s+(r2),r2-ri) < r, (Case I) 
/i2(s-(r2),s+(r2),r2;s_(r2)-s_(ri),s+(r2)-S-H(ri),r2-ri) < r, (Case II) 
where /ij’s are the functions in Theorem | 5 . 5 t and 

T := min{( 5 , 5(s+(r2) - s_(r2))/4}, (Case I) 
r := min{ 5 , (5(s 2 — ■Si)/ 4 }. (Case II) 

Now, fix any ri G and let B be the n x n matrix determined 

through Steps 2 and 3 in terms of any given {po,f3o) £ S = Sr^^r2- Let 
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= po + toqo and p- = po + soqo from Step 2 ; then p± and A(p±) fulfill 
the conditions in Theorem 5.5 So this theorem implies that there exists a 
vector G such that 

max{|p° -p_I,Ip^ -p+l,IA(p°) -A(p_)l, |^(p+) - ^(p+)|} < r, 
where p\ = kC,^, p^_ = IC,^, and A{p\) = r2C°- Using ( 5 . 7 ) and ( 5 . 8 ), 
\p+-K^\< 5 , \p.-T(f\< 5 , 

|gg _ ^0| ^ | P+ -P- _ ^0| < l(p+ -P-) - (fc- 0 C°l + 1(^-0 - (to - So)| 


to — So 

^ 2t + ||p° I - |p+ -p_|| ^ 
~ to — So 


to — So 
4t 

to — So 
.0 


< 5 , 


I I l^(p+)-^(p-)l ^ l^(p+)-^(p+)l + l^(p-)-^(p-)l ^ ^ 

l 7 o| = I- - - 1 < - - -< 0. 

to — So to — So 

Since det(i?) = DET(p+,p_, g'Oj 7 o) and |DET(fcC‘^, 0 )| > d, it follows 

from ( 5 . 10 ) that 

I det( 71 )| > d/2 > 0 . 

The proof is now complete. □ 


5 . 2 . Relaxation of Vuj{z) G Kq. The following result is important for the 
convex integration with linear constraint; the function p determined here 
plays a similar role as the tile function g used in | 43 l 03 ] . For a more general 
case, see | 37 l Lemma 2 . 1 ]. 

Lemma 5 . 7 . Let Ai, A2 > 0 and pi = —Airy, 772 = X2'n with 

"“(h®.! 7 )’ M = 1 . 7 -<I = 0,M0, 

Let G C be a bounded domain. Then for each e > 0, there exists a 

function u = G with supp(a;) CC G that satisfies 

the following properties: 

(a) div?/; = 0 in G, 

(b) |{ 2 ; G G I Vlo{z) i {jyi, 772}}! < e, 

(c) dist(Va;(2;), [771,772]) < e for all z £ G, 

(d) ||a;||£,oo(-G) < e, 

(e) (p{x,t) dx = 0 for all t G'R. 

Proof. The proof follows a simplified version of | 37 [ Lemma 2 . 1 ]. 

1 . We define a map V: G^(M"''''^) —)■ G'^(M"''''^;by setting V{h) = 
{u,v), where, for h{x,t) G G^(M’^+^), 

u(x, t) = q ■ Dh{x, t), v{x, t) = -{'y'S> q — q 'S) 'y)Dh(x, t). 

b 

We easily see that V{h) = {u,v) G G“(M"'+^; supp(T’(/i)) C supp(/7), 

divu = 0 , and u{x,t) dx = 0 for all t G M, for all h G G“(M”'''^). For 
h{x,t) = f{q ■ X + bt) with / G G°°{R), w = {u,v) = V{h) is given by 
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u{x,t) = f'{q ■ X + bt) and v{x,t) = f'{q ■ x + bt)^, and hence Vw{x,t) = 
f"{q ■ X + bt)r]. We also note that V{gh) = gV{h) + hV{g) and hence 

( 5 . 11 ) Vr{gh) = gVV{h) + hVV{g) + B{Vg, Vh) \/ g, h e 

where B{S/g,'Vh) is a bilinear map oiS/g and Vh] so \B{S/h,'Vg)\ < ClV/il | Vg'| 
for some constant C > 0 . 

2. Let Ge CC G be a smooth sub-domain such that |G \ Ge| < e/2, and 
let pe G G“(G) be a cut-off function satisfying 0 < pe < 1 in G, pe = 1 on 
Ge- As G is bounded, G C {(x,t) \ k < q ■ x + bt < 1 } ior some numbers 
k < 1 . For each r > 0 , we can find a function f^- G C^{k, 1 ) satisfying 

-Al < f” < A2, \{s G {k,l) I //(s) ^ {-Ai, A2}}| < r, ||/r||L°° + ||/rl|L“ < T. 

3 . Define io = (<p, V’) = 'P{pe{x,t)hr{x,t)), where hr{x,t) = friq ■ x -|- 
bt). Then \\hr\\c^ < C\\fr\\c^ < Ct, uj G G“(M’^+^;supp(a;) C 
supp(pe) CC G, and (a) and (e) are satisfied. Note that 

|w| < \pe\\V{hr)\ + \hr\\V{pe)\ < G^T, 

where Ge > 0 is a constant depending on ||pe||ci(G)- So we can choose a 
Ti > 0 so small that (d) is satisfied for all 0 < r < ri. Note also that 

{z G G I Vu:{z) ^ {pi,P2}} C (G\Ge)U{z G G^ | f”{q-x + bt) ^ {-Ai, A2}}. 

Since \{z G G, | f"{q-x + bt) ^ {-Ai, A2}| < N\{s G {k,l) \ f”{s) ^ 
{—Al, A2}}| for some constant > 0 depending only on set G, there exists 
a r2 > 0 such that 

|{z G G| Va;(z) ^ {pi,P2}}| <^+Nt <e 
for all 0 < r < T2. Therefore, (b) is satished. Finally, note that 
p^W{K{x,t)) = Pef”{q ■ X -h bt)g G [pi,P 2 ] in G 
and, by d^ , for all z = {x, t) G G, 

\^i0{z) - p,VV{hr{x,t))\ < \K\\VV{p,)\ + \B{Vhr,Vp,)\ < C',T < e 

for all 0 < r < Ts, where G' > 0 is a constant depending on ||pe||c 2 (G)) 
and Ts > 0 is another constant. Hence (c) is satisfied. Taking 0 < r < 
min{Ti, r2, Ts}, the proof is complete. □ 

We now state the relaxation theorem for homogeneous differential inclu¬ 
sion Vu){z) G Kq in a form that is more suitable for later use; we restrict 
the inclusion only to the diagonal components. 

Theorem 5 . 8 . Let 0 < r2 < cr(so) (Case I), 17(52) < r2 < (x{si) (Case 
II), and let I2 = ^ G ( 0 ,r 2 ) (Case I), I2 = lr2 G (<7(52), 72) (Case II) be 
some number, depending on r2, from Theorem \ 5 . 6 [ Let I2 < ri < r2, and 
let 1 C be a compact subset of S = Sr-,^,r2- Q ^ I be a box in Then, 

given any e > 0 , there exists a 6 > 0 such that for each box Q x L C Q x I, 
point {p, (d) G 1 C, and number p > 0 sufficiently small, there exists a function 
u = ((p, i/) G G“(Q X I;]R^+"') satisfying the following properties: 
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(a) divi/) = 0 in Q X I, 

(b) (p'+D(p{z), l3'+^l>tiz)) G S for all z € Qxl and \{p', /3')-{p, 13)\ < 6, 

(c) \\ui\\l^{qxI) < Pi 

W IqxI\ f^ + “ ^(P + Dip{z))\dz < e\Q x I\/\Q x /|, 

IqxI P + Mz)), ~^)dz < e\Q X I\/\Q x /|, 

(f) Jq p(x, t)dx = 0 for all t £ I, 

(9) \\Pt\\L°°{QxI) < Pi 

where A = Ar^^r 2 C is the set defined by 

= WiA{p'))\\p\ e [s-(ri),s_(r2)] U [s+(r2),s+(ri)]}, (Case I) 


■^= {{PiA{p'))\W\ e ['S-(ri),S-(r2)] U [s+(ri),s+(r2)]}. (Case II) 

Proof. By Theorem |5.6[ there exist finitely many open balls Bi, - ■ ■ ,Bn CC 
S covering K, and functions qi : Bi ^ 7* : —>■ M”", ti^± : —)• M 

(1 < i < N) with 7i • gi = 0 and U- < 0 < on Bi such that for each 
p c 


e = 


B /? 


G R{Fq) with {p,l3) £ Bi, we have 

^ T ti^^pi £ Fzizi 
QiiPi P) 


where = (.,±(p,/3), Vi = (/,) “d M 0 ie 

arbitrary. 

Let 1 < z < iV. We write f,i = ii{p,fd) = (^q ^ R(Fo) for (p, ft) £ 

Bi C S, where O is the n x n zero matrix. We omit the dependence on 
{p, (3) £ Bi in the following whenever it is clear from the context. Given any 
p > 0, we choose a constant bi with 

0 < < min- - -. 

Bi li,+ I'i,— 

With this choice of b = bi, let pi be defined on Bi as above. Then 




Pi,± Cj,± 
Bi,± /3i,± 


■— f,i + ti,±Vi £ -^±1 


= Aiei,+ + (l-\)4-, Xi = - -^5^G(0,1) on 

^i,+ ti^— 

By the definition of R{Fq), on Bi, both + (1 — T)^i- and 

= (1 — r)^j,+ + T^j,- belong to R{Fo) for all r G (0,1). Let 0 < 
r < mini<j<Ar ming. min{Aj, 1 — Aj} < ^ be a small number to be selected 
later. Let A' = ^ on Bi. Then A' G (0,1), = A'4^+ + (1 - 

on Bi. Moreover, on Bi, = (1 — 2T)(^j_+ — f,i-) is rank-one, 

[f.l_Al+] C (^i_,^i,+) C R{h), and 

CT < - ^7+1 = - ^i,-| = r|^i,+ - ^i-\ = r(ti,+ - ti-)\pi\ < Ct, 
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where C = maxi<j< 7 vmaxg^. 

= c> 0. By continuity, Hr = U(p,/3 )eB,,i<i<Af[^J-(f'’^)] ^ 

pact subset of R{Fq), where R{Fq) is open in the space 


Bn = 


p c 

B (3 


HB = 0 


by Lemma 5.3 and Theorem |5.6[ So dr = dist(LfT-, 9|so.R(-^o)) > 0, where 
9|so is the relative boundary in Eq. 

Let = -Xippi =-A'(l - 2r)(L,+-L _)r/i, pi ^2 = Xi^ 2 Pi = (1-A')(1- 
2r)(L,+ - ti-)pi on Bi, where Aj,i = r(-ti,+) + (1 - T){-ti-) > 0, Ai ,2 = 
(1 — r)L^+ + Tti- > 0 on Bi, and r > 0 is so small that 

min minAoi.>0 (fe = 1, 2). 
i<j<N Bi 


Applying Lemma 5.7 to matrices pi^i = pip{p,f3), pi ^2 = 'ni, 2 {p, P) with a 
fixed {p, (3) G Bi and a given box G = Q x I, we obtain that for each 
p > 0, there exist a function ui = {ip, 'll;) G C^{Q x /; and an open set 

Gp CC Q X I satisfying the following conditions: 

(5.12) 

(1) divi/i = 0 in Q X /, 

( 2 ) |(Q X I) \ Gp\ <p- ^i + Vuj{z) G {il_, Cl+} for all z G Gp, 

(3) + Vlo{z) G [4"-, e*"+]p for all z G g X I, 

(4) ||a;||2,tx>(Qx/) < Pi 

(5) Jg p{x, t) dx = 0 for all t £ I, 

.(6) ||<Pt||L°°(Qx/) < 2p, 

where [^l_,^l^]p denotes the p-neighborhood of closed line segment 
Here, from (5.12 3), condition (5.12 6) follows as 

\ipt\ < \ci,+ - Ci-\ + p = (L,-h - ti-)\bi\ + p <2p in g X /. 


Note (a), (c), (f), and (g) follow from (5.12), where 2p in (5.12 6) can 
be adjusted to p as in (g). By the uniform continuity of A on the set 
J = {p' G M” I Ip'I < s+(ri) (Case I), \p'\ < s+(r 2 ) (Case II)}, we can 


find a (5' > 0 such that \A{p') — A{p'')\ < 


3|Qx/| 


whenever p', p" G J and 


\p' — p''\ < 5'. We then choose a r > 0 so small that 

Gt < 5', G\Q X i\T < |. 


Next, we choose a <5 > 0 such that <5 < ^. If 0 < p < <5, then by (5.12 1) 
and (5.12 3), for all z G g x / and \{p',(3') — {p,f3)\ < 5, 

^i{p, /30 + Vw(z) G Eo, dist(^i(p', (3') + Vuj{z),Hr) < dr, 


and so ^i{p',f3') + Vuj{z) G R{Fo), that is, {p' + Dip{z),l3' + i>t{z)) G S. 
Thus (b) holds for all 0 < p < <5. In particular, {p + Dip{z),j3 + '4’t{z)) G S 
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and so \p + D(p{z)\ < s+(ri) (Case I) 


1/3 + '0t(z)| <r 2 for all z G Q X I, by (i) of Theorem 5.6 


\p + D(p{z)\ < s+{r 2 ) (Case II) and 
Thus 


'Qxl 


1/3 + i/Zt - A{p + D(p)\dz 


< [ \/3 +'ijjt - A{p +D(p)\dz + {r 2 +M^)p 

JGp 

< IQ X /| max{|/3/)± - A{pl^)\} + (r 2 + M„)p 

< CIQ X /|r + IQ X /| max{|^(pi,±) - A(p^)|} + (r2 + Mu)p 


< 


2e|Q X I\ 
3|Q X i\ 


+ (r2 + M„)p, 


where = <7('So) (Case I), = cr(s 2 ) (Case 

’ \^h± Pi,±J 

II). Thus, (d) holds for all p > 0 satisfying (r 2 + Mfj)p < 5 !^^. Similarly, 

o|(c,/X i I 


iQxI 


dist((p + Dp{z), 13 + 'iptiz)),A)dz 


< 


/ max|(p[±,/3,’(±) - {pi,±, I3i,±)\dz + 2 (r 2 + M^)p 

JGp 


< C\Q X I\T + 2{r2 + M„)p 
3|Q X I\ 


therefore, (e) also holds for all such /o > 0 for (d). 

We have verified (a) - (g) for any {p, j3) G Bi and 1 < i < A^, where 5 > 0 
is independent of the index i. Since Bi, - ■ ■ ,Bn cover /C, the proof is now 
complete. □ 


6. Boundary function and the admissible set U 


Assume hi and uq satisfy (1.6) and (2.1). 


6.1. Boundary function <h. We first construct a suitable boundary func¬ 
tion = (u*,v*) in each case of the profile (t(s). Recall that the goals of 


Cases I and II are to prove Theorems 3.1 and |3.5| under Hypotheses (PM) 
and (H), respectively. 


Case I. In this case, H is assumed to be convex. Let 0 < r = r 2 < (t(Mo), 
and let / = Q 2 ^ (0) ^) be some number determined by Theorem 5.6, Choose 
any f = ri G {I, r). 
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Case II. In this case, we assume |Dno(xo)| £ (■Si,S 2 ) some xq G VL. 
Let ct ( s 2 ) < r = r 2 < a{MQ), and let I = lr 2 G (o'(s 2 ),r) be some number 
determined by Theorem 5.6 Choose any f = ri G (l,r). 


We now apply Lemma 4.2 (Case I), Lemma 4.3 (Case II) to deter¬ 
mine functions ct, / G (^^([OjOo)) (Case I), d, / G (^^^“([O, oo)) (Case II) 
satisfying its conclusion. Also, let A(p) = f{\p\‘^)p {p G M”). Then 

Lemma 6 . 1 . We have 

(p,A(p))gS V s-(ri) < IpI < s+(r 2 ), 
where S = Sri,r 2 the set in Lemma \5.3[ 

Proof. Let s = \p\, r = cj(s) and C = p/ \p\, so that s_(ri) < s < s+(r 2 ), 


C G ^ and A(p) = rf. By Lemma 
s < S-i-(r) and ri < r < 

A{p±) = rC, = I5±. Define ^ = 


4.2 


4.3 


(Case II), 


S-{r) < s < S-i-(r) and ri < r < r 2 . Set p± = s±{r)C, and f3± = rC^. Then 

O A(p)) ““ = (o L) ■ 

^ = A^+-t-(l—A)^_ for some 0 < A < 1. Since G F± and rank(^+—^_) = 1, 
it follows from the definition of R{Fo) = i2(Frj^r2(0)) that ^ G (^_,^+) C 
R{Fo). Thus, by Lemma [s^ {p,A{p)) G 5. □ 

By Lemma fd . 2| (Case I), Lemma [4.3| (Case II), equation ut = div(A(iAri)) 
is uniformly parabolic. So by Theorem |4.1t the initial-Neumann boundary 
value problem 


( 6 . 1 ) 


= div(A(iAu*)) in D't 
du*/dn = 0 on dPt x (0,r) 

tt*(x,0) = uo{x), X G D 


admits a unique classical solution u* G moreover, only in 

Case I, it satisfies 

\Du*{x,t)\ < Mq V(x,t) G Qt- 


From conditions (1.6) and (2.1), we can find a function h G C^^“(D) 
satisfying 

Ah = uq in Q, dh/dn = 0 on dQ. 

Let uo = Dh £ C'^'''“(D; M”) and define, for {x,t) G Dr, 

(6.2) V* {x,t) = vo{x) + f A{Du*{x, s)) ds. 

Jo 

Then it is easily seen that <I> := {u*,v*) G (7^(0^;satisfies (2.2); that 
is, 


(6.3) 


u*{x, 0) = uo{x) (x G D), 
divu* = u* a.e. in Dr, 

■ n|an = 0 Vt G [0,r], 
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and so is a boundary function for the initial datum no- 
Next, let 

^={{P:Mp)) I \P\ e [0>s-(n)]}, (Case I) 

^ = {{p,Mp)) I \P\ e [ 0 ,s-(n)] U [s+(r2),max{Mo,S2}]} . (Case II) 
Then we have the following: 


Lemma 6 . 2 . 


{Du*{x, t),vl{x, t)) £ S U T V (x, t) £ Qt- 


Proof. Let {x,t) £ and p = Du*{x,t). 

Case I: In this case, \p\ < Mq < s+(r 2 ). If \p\ < s_(ri), then A{p) = A{p) 
and hence by ( |6.2| ) 

{Dv*{xA),vl{x,t)) = {p,A{p)) = {p,A{p)) £F. 


If S-(ri) < \p\ < Mq, then by Lemma 6.1 and (6.2) 


{Du*{x,t),Vtix,t)) = {p,Mp)) e •S- 


Case II: If \p\ < s_(ri) or s+(r 2 ) < \p\ < Mq, then A(p) = A{p) and 
hence as above 


{Du*{x,t),vt{x,t)) = ip,A{p)) = ip,A{p)) £T. 

If S-(ri) < \p\ < s+(r 2 ), then as above 

{Du*{x,t),vl{x,t)) = {p,Mp)) e S. 

Therefore, for both cases, {Du*,v^) £ S L) P' in il,T- Cl 

Definition 6 . 3 . We say a function u is piecewise in Q.t and write u £ 
Cpiece(^r) if there exists a sequence of disjoint open sets in Q.t 

snch that 

u £ C^{Gj) Vj G N, |flT \ ^f=iGj\ = 0. 

It is then necessary to have \dGj\ = 0 Vj G N. 


6 . 2 . Selection of interface. To separate the space-time domain fly into 
the classical and micro-oscillatory parts for Lipschitz solutions, we assume 
the following. 


Case I. Observe that 


|{(x, t) £ Qt\ \Du*{x, t)| = s_(r)}| >0 
for at most countably many f G (0, f). We fix any f G (0, f) with 
|{(x,t) G Qt I \Du*{x,t)\ = s_(f)}| = 0, 

and let 

= {(®) t) G O-r I \Du*{x, t)| < s_(f)}, 

= {{x, t) £ I \Du*{x, t)| > s_(f)}, 
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SO that and are disjoint open subsets of Qt whose union has full 
measure |n'r|. Define 

Dg = {(x, 0) I X G D, |L)uo(x)| < s_(f)}; 

then Dg C 

Case II. Observe that 

|{(x,t) G Dr I |i:>u*(x,t)| = s_(fj)}| >0 (j = 1,3) 

for at most countably many fi G (< 7 ( 52 )) f), ^3 G (r, cj(si)). We fix any two 
fi G {a{s 2 ),r), fs G (r,cj(si)) such that 

|{(x,t) G Drl |i:>u*(x,t)| = s-(fj)}| = 0, (j = 1,3) 

and let 

Dr = {(x, t) G Dr | \Du*{x, t)| < s_(fi)}, 

Df. = {(x,t) G Dr I \Du*{x,t)\ > s+(f 3 )}, 

and 

Dy = {(x,t) G Dr I s-(ri) < \Du*{x,t)\ < s+(f 3 )}, 
so that these are disjoint open subsets of Dr whose union has full measure 
|Dr|. Define 

Dgi = {(x,0) I X G D, |Duo(x)| < s_(fi)}, 

Dq 3 = {(x,0) I X G D, |Iluo(x)| > s+(f 3 )}; 
then Dg-’ C 9D^ (j = 1,3). 

6 . 3 . The admissible set U. Let m = ||rti ||roo(Dy) + 1 - We finally define 
the admissible set U and approximating sets (e > 0) as follows. Recall 
0<f<f = ri<r = r 2 < ct(so) (Case I), cr(s 2 ) < n < t < r = r 2 < 
f3 < ct(si) (Case II). 

Case I. 

U = n W^i“(Dr) I u = u* in D;^, ||ut||roo(Q,j,) < m, 

3v G n Wji°°(Dr;M'^) such that u = u* in D^, 

divu = u and {Du, vt) G S U D a.e. in Dr, and 

{Du, vt) G S U B a.e. in D|,|, 

Ue = ^uGh(\3vG Cp^^^^ n VLj*°°(Dr; M”) such that u = u* in D^, 
div V = u and {Du, vt) G S U D a.e. in Dr, 

{Du,vt) gSuB a.e. in D^, \vt — A{Du)\dxdt < e|Dr|, and 

f ^2 dist{{Du,vt),C)dxdt < €|D^||, 

where 

B = Br,f,r = {{p,A{p)) I IpI g [s-{f),s-{r)]}, 

C = Cf,f,r = {{p,A{p)) I \p\ G [s-{r),s-{r)] U [s+(r), s+(f)]}. 
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Case II. 

U = Cpjgce n Wlf°{Q.T) I -U = u* in U VL^, \\ut\\(^U t) < rn, 

Bug n Wlf°{^lT] M”) such that u = u* in U 

div V = u and {Du, vt) G 5 U a.e. in VLt, and 

{Du, Vt) £ S \J B a.e. in 


Z^e = |tt G I 3 u G n Wlf^{DT', M”") such that u = u* in Qp U Qp, 

divu = u and {Du, vt) G S U D a.e. in Qt, 

{Du,vt) gSVJB a.e. in \vt — A{Du)\dxdt < 61^7^1) and 

f ^2 dist{{Du,vt),C)dxdt < e|r2|.||, 

where 

B = Bf^,r,r,f 3 = {{p,A{p)} I \p\ G [s_(fi),s_(f)] U [s+(r),s+(f 3 )]}, 


C = Cfi,f,r,f 3 = {{p,Mp)) I \p\ e [s_(ri),s_(r)] U [s+(f), s+(r 3 )]}. 

Observe that for both cases, AU B = C and B <Z T, where A is as in 
Theorem |5.8[ Also for both cases, as in the proof of Lemma 6.2, it is easy 
to check that {Du* ,v^ ) G 5 U .6 in 

Note that one more requirement is imposed on the elements of in both 
cases than in the general density approach in Subsection 2.2 As we will 


see later, such a smallness condition on the distance integral is designed to 
extract the micro-structured ramping of Lipschitz solutions with alternate 
gradients whose magnitudes lie in two disjoint (possibly very small) compact 
intervals; this occurs only in and the solutions are classical elsewhere. 

Remark 6 . 4 . Summarizing the above, we have constructed a boundary 
function $ = {u*,v*) G {CIt',^^^^) for the initial datum uq with u* G Id; 
so U is nonempty. Also U is a bounded subset of Wlf°{QT)^ since 5 U T” is 
bounded and ||itt ||loo(Oj,) < m for all u GlA. Moreover, by (i) of Theorem 


5.6 and the definition of T, for each u GlA, its corresponding vector function 
V satisfies \\vt\\L^(nj,) <r = r2 (Case I), ||ut< max{r = r 2 ,a{Mo)} 
(Case II) ; this bound in each case plays the role of a fixed number R > 
0 in the general density approach in Subsection |2.2[ Finally, note that 
s_(ri) < \Du*\ < s+{r2) on some nonempty open subset of Qt and so 
A{Du*) A A{Du*) on a subset of fl-r with positive measure; hence u* itself 


is not a Lipschitz solution to (1.2). 


In view of the general existence theorem. Theorem 2.2 it only remains 


to prove the L°°-density of lA^ in U towards the existence of infinitely many 
Lipschitz solutions for both cases; this core subject is carried out in the next 
section. 
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7. Density of lie in 7/: 

Final step for the proofs of Theorems 13.11 and 13.51 

In this section, we follow Section to complete the proofs of Theorems 
13.11 and 13.51 


7.1. The density property. The density theorem below is the last prepa¬ 
ration for both cases. 


Theorem 7.1. For each €> 0,lle is dense inU under the L°°-norm. 

Proof. Let u G U, rj > 0. The goal is to construct a function u G He such 
that \\u — < r]. For clarity, we divide the proof into several steps. 

1. Note that u = u* in (Case I), in (Case II), ||L°o(Hy) < 

m — To for some tq > 0, and there exists a vector function v G n 

lTj*°°(D'r; M"') such that u = u* in (Case I), in U (Case II), 
divu = u and {Du,vt) G S U T a.e. in Dy, and {Du,vt) G S U B a.e. in 
Dy. Since both u and v are piecewise in Dy, there exists a sequence of 
disjoint open sets {Gj]Y=i ™ with \dGj\ = 0 such that 

u G G^{Gj), V G C^(Gj;M") Vj > 1, |D| \ = 0. 

2. Let j G N be fixed. Note that {Du{z),vt{z)) G SVJB for all z = (x,t) G 
Gj and that Hj := {z G Gj \ {Du{z), vt{z)) G dS} is a (relatively) closed 
set in Gj with measure zero. So Gj := Gj \ Hj is an open subset of Gj with 
l^il = \Gj\, and {Du{z), vt{z)) G S Gi B ioi nW z G Gj. 

3. For each r > 0, let 

Qr = {(p,/3) G S I dist((p,/3),d5) > r, dist((p,/3),C) > r} ; 

then 

s = {Gr>oQr) u {(p,/3) G S I dist((p,/3),C) = 0}. 

Since B C C, we have 

I \vt{z) - A{Du{z))\ dz 
Jg, 

= lim / \vt{z) — A{Du{z))\dz, 

I (Du(z),vt(z))£g.r} 


iGi 


dist{{Du{z), vt{z)), C) dz 


= lim / 

^^0+ .J{zeGj I (Duiz),vt(z))&GT} 


dist((Du(z), vt{z)),C) dz] 


thus we can find a > 0 such that \dOj\ = 0, 
(7.1) \vt{z) - A{Du{z))\dz < 
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and 

(7.2) 


dist{{Du{z), vt{z)), C) dz < 


3 • 


I^tIj 


where Fj = {z € Gj \ {Du{z), vt{z)) ^ Qr^} and Oj = Gj \ Fj is open. Let 
J be the set of all indices j G N with Oj / 0. Then for j 0 J, Fj = Gj. 

4. We now fix a j G J. Note that Oj = {z G Gj \ {Du{z), vt{z)) G Qt^} 
and that JCj := Qr^ is a compact subset of S. Let Q C M"" be a box with n C 
Q and I = (0, T). Applying Theorem 


5.8 


to box Q X I, fCj CC S = S. 


and e' = we obtain a constant 5 j > 0 that satisfies the conclusion of 

the theorem. By the uniform continuity of A on compact subsets of M”, we 
can find a, 9 = > 0 such that 

(7.3) \A{p) - A{p')\ < ^ 

whenever |p|, \p'\ < s and \p — p'\ < 0 , where the number s > 0 will be 
chosen later. Also by the uniform continuity oi u, v and their gradients on 
Gj, there exists a Uj > 0 such that 

\u{z) — u{z')\ + |Vu(z) — Vu(z')| + \v{z) — v{z')\ 


(7.4) 


+|Vu(z) - Vv{z')\ < min{^, ^,9} 


whenever z,z' G Gj and |z — z'| < Vj. We now cover Oj (up to measure 
zero) by a sequence of disjoint boxes {Q*- x in Oj with center z*- and 


diameter < Uj. 


5. Fix an i G N and write w = {u,v), ^ = 


p 

B 


Du{zi] 


Mz}) 


Dv{z)) vt{z] 


= Vw(zj) = 
5.8 


since 


By the choice of 6j > 0 in Step 4 via Theorem 

jj -rv-jV - 

Q^j X I J C Q X I and {p, j3) G K,j, for all sufficiently small p > 0, there exists 
a function w) = (</?), V’!) G 0'^{Q^, x /G]R^+"') satisfying 


(a) 

(b) 

(c) 

(d) 

(e) 

(f) 

(g) 


div?/;*- = 0 in Q* x /*, 

{p' + Dip^-{z),l5' + {'iljj)tiz)) G S for all z € x F- 
and all |(p',/3') - (p,/3)| < 5j, 

Jq^xT 1/3 + i^M^) - Mp + D^iiz))\dz < d\Q) x^\/\Qxi\, 
/Q5x7jdist((p + 7)(/?j(^;),^ + {'il^i)t{z)),A)dz < e'lQ) x /j|/|Q x I|, 
fgi ip^j{x, t)dx = 0 for all t £ Ij, 

\\iP))t\\L°°{Q)xr.) < P^ 

C is as in Theorem 


'-J 3 

where the set A = A. 

5 . 


ri,r2 


5.8 


Here, we let 


4.4 


0 < p < min{To, 2 ^, where > 0 is the constant in Theorem 

and O is the product of On and the sum of the lengths of all sides of Q. 
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From xD) = (f), we can apply Theorem 

to obtain a function = TZip'j G C^iQ^j x Ij; M") n 
that div g't = (p\ in Q*. x P- and 


4.4 to (f'j on Q'j X /] 
’“(Q* X T;M") such 


(’^•5) ll(5j)t|li,°°(Qix/j) — ^\\i^])t\\L°°(QixPj) — (by (g)) 

6. As \vt — A{Du)\, dist((T>u, ut), C) G L°°{Qt), we can select a finite 
index set X C J x N such that 


(7.6) 


'U(. 




\vt{z) - A{Du{z))\dz < 

o 


(7.7) [ dist{{Du{z),vt{z)),C)dz < 

We finally define 

{u, v) = {u, v)+ ™ 

Old el 

As a side remark, note here that only finitely many functions + gj) 

are disjointly patched to the original {u, v) to obtain a new function (u, v) 
towards the goal of the proof. The reason for using only finitely many pieces 
of gluing is due to the lack of control over the spatial gradients + 9j), 

and overcoming this difficulty is at the heart of this paper. 

7. Let us finally check that u together with v indeed gives the desired 
result. By construction, it is clear that u G n v G Clbggg n 

IK"') and that u = u = u* and v = v = v* m (Case I), in 
(Case II). By the choice of p in (g) as p < tq, we have \\ut\\(Q j.) < 
m. Next, let {j,i) G X, and observe that for z & x I), with {p, fi) = 
{Du{Zj),vt{Zj)) G Gtj, since \z — z'j\ < Ij < uj, it follows from (7.4) and (7.5) 
that 

\{Du{z),vt{z) + {gj)t{z)) - < Sj, 

and so {Du{z),vt{z)) G S from (b) above. From (a) and divg^ = for 
zeQ)x /*, 

divi}(z) = div(u + fij + g]){z) = u{z) + 0 + fi^j{z) = u(z). 
Therefore, u G U. Next, observe 

\vt — A{Du)\dz 

\vt — A{Du)\dz = / \vt — A{Du)\dz 
\vt - A{Du)\dz + f 

iii Jui 


/ \vt — A{Du)\dz = / \vl — A{Du*)\dz + 






\vt — A{Du)\dz 


''~^(j,i)e(Jxn)\xQ'jXl 
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— • -^1 + -^2 + -^3) 


/n2 


disi{{Du,vt),C)dz = / dist{{Du,vt),C)dz 


+ 




dist{{Du,vt),C)dz + / dist{{Du, vt),C)dz 


= :l! + ll + ll 


From (7.1), (7.2), (7.6) and (7.7), we have if + /| < y (A; = 1, 2). Note 
that for {j,i) G I and z £ x Ij, from (7.4), (7.5) and (g), 

\vt{z) - A{Du{z))\ = \vt{z) + {i^^Mz) + {g"Mz) - A{Du{z) + D^p){z))\ 


< \vt{z) - Vt{z^j)\ + \vt{z^j) + - A{Du{z^j) + Dlp){z))\ 

+\{9''j)t{z)\ + \A{Du{z^j) + D^p){z)) - A{Du{z) + Dip){z))\ 
< g + \vt{z^j) + {^^t{z) - A{Du{z)) + Dip){z))\ 

+ \A{Du{z^j) + Dip{j{z)) - A{Du{z) + DLp'^j{z))\. 
Similarly, since ^ C C, we have 


disi{{Du{z),vt{z)),C) 

< I + dist((Zl'«( 2 :*) + Dip]{z),vt{z''j) + {ip))t{z)),C) 

< I + dist{{Du{Zj) + Dip{j{z),vt{zj) + {^pj)t{z)), A). 


From (b) and (i) of Theorem 5.6, we have \Du{Zj) + D(p{j{z)\ < s+(ri) =: s 
(Case I), \Du{zj) + Dip{j{z)\ < s+(r 2 ) =: s (Case II). As {Du{z),vt{z)) £ 
S, we also have \Du{z) + Zl(/?t(z)| = \Du{z)\ < s, and by (7.4), \Du{zj) — 
Du{z)\ < 9. From (7.3), we thus have 

\A{Du{z)) + Dif]{z)) - A{Du{z) + D^]{z))\ < 


Integrating the two inequalities above over Q*- x Ij, we now obtain from (d) 
and (e), respectively, that 



el^Tl \Q) X ^ 
12 IQ X i\ 




d\si{{Du{z),vt{z)),C)dz < ^|Q* x/j| + 


elOrl IQj X P- 
12 IQ X i\ 


<\\Q 


• X /'I 


thus < ||f2j’|, and so + if + -^3 — cI^tI) where k = 1,2. Therefore, 
u £ lAe- Lastly, from (c) with p < rj and the definition of u, we have \\u — 
< '9- 

The proof is now complete. □ 
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7.2. Completion of the proofs of Theorems |3.1| and |3.5| . Unless specif¬ 
ically distinguished, the proof below is common for both Case I: Theorem 
13.11 and Case II: Theorem [ 


Proofs of Theorems \3.1\ and 3.5, We return to Section As outlined in 
Remark 6.4 Theorem |7.1| and Theorem 2.2 together give infinitely many 


Lipschitz solutions u to problem (1.2) 


We now follow the proof of Theorem 2.2 for detailed information on such 


a Lipschitz solution u £ Q to (1.2). Here Du is the a.e.-pointwise limit 


of some sequence Duj, where the sequence uj G Ui/j converges to u in 
Since Uj = u* in (Case I), in (Case II), we also have 

u = u* G (Case l),u = u* G C2+“d+«/2(f^^ y ) (Case 

II) so that 

ut = div{A{Du)) and \Du\ < s_(f) in (Case I) 

ut = div(A(i7u)), \Du\ < s_(fi) in 0^, 

and I Du I > s+(f 3 ) in (Case II) 

Note {vj)t vt in L^(n'r;M"'), where Vj is the corresponding vector func¬ 
tion to Uj and v G VU^’^((0, T); 7.^(0;M"-)). From (2.8), we can even deduce 
that {vj)t —)• Vt pointwise a.e. in Qt- On the other hand, from the definition 


oiU 


i/i’ 




dist{{Duj, {vj)t),C) dxdt < —)• 0 as j —>■ oo; 


thus {Du,vt) G C a.e. in yielding |5| -|- |L| = 

For the remaining assertions, we separate the proof for each case. 

Case I. If |L| = 0, then |Du| < s_(r) a.e. in H-r; so u is a Lipschitz 


solution to (6.1) with monotone flux A{p). Thus, by Proposition 

This contradicts the fact that ||Du 


have u = u* in fir 
11 Duo 11 LOO = Mo > S_(r). Thus |L| > 0. 

Case II. Suppose |L| = 0. Then 


•,o)|| 


2.3 


we 


L°°(n) 


|Du| G [0, s_(r)] U [s+(r 3 ),oo) a.e. in 

Now, modify the profile cr(s) so as to obtain a function a G C'^+“([0, oo)) 
satisfying 

a(s) = o-(s),_ s G [0,s_(r)] U [s+(f 3 ),oo), 

6 < ^'{s) <0, 0 < s < oo, 

for some constants 0 > (9 > 0, and let /(s) = a^{^/s)/^/s (s > 0), /(O) = 
/(O), and A(p) = /(|pp)p (p G MA). Then the functions u G Q are all 
Lipschitz solutions to the problem of type (1.2) with monotone flux A(p), a 


contradiction to the uniqueness by Proposition 2.3 Thus |L| > 0. 
Next, suppose |5| = 0. Then 

|Du| G [0, s_(fi)] U [s+(r),oo) a.e. in rir- 
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So we get a contradiction similarly as above by obtaining a function a G 
(27i+o([o,oo)) satisfying 

a(s) = cr(s),_^ s e [0,s_(fi)] U [s+(f),oo), 

6 < a'{s) <0, 0 < s < oo, 


for some constants 0 > 0 > 0. Thus |5| > 0. 


□ 


7.3. Proof of Corollary | 3.6 Recall that this corollary is under Case II: 
Hypothesis (H). Let uq G with Duq ■ n|aQ = 0. The existence 

of infinitely many Lipschitz solutions to problem (1.2) when |Z)tto(xo)| G 


( 5 ^, 52 ) for some xq G H is simply the result of Theorem 3.5 So we cover 
the other possibilities here. 

Assume ||L)uo 11^00(^2) < s^. Fix any two numbers cr(s 2 ) < ri < r 2 < 
and let d, / G C'^^“([0, 00 )) be s ome functions from Lemma 
Using the flux A{p) = .f(|pP )iP, Theorem 
C 2 +“’i+«/ 2 (j 7 ^) to problem (Q. If \Du* 
si in Qj', then u* itself is a Lipschitz solution to (1.2). 


4.1 


4.3 


gives a unique solution u* G 


stays on or below the threshold 
Otherwise, set 


_ ^l+g-(n) 


s = 


and choose a point (x, t) G VLt such that 
\Du*\ <s in Q X (0,t), \Du*{x,F)\ G 

Regarding ui(-) := G satisfying Dui • n|af 2 = 0, as a new 


initial datum, it follows from Theorem 3.5 that problem (1.2), with the 


initial datum ui at time t = t, admits infinitely many Lipschitz solutions u 
in H X {t,T). Then the patched functions u = Xnx( o.F)^ * + Xnx[t,T)U in 
become Lipschitz solutions to the original problem ( |1.2[ ). 

Lastly, assume min q \Du o\ > Let ri, r 2 and A{p) be as above, and let 
u* be the solution to (4.3) corresponding to this flux A{p). If \Du*\ stays 


on o r above the threshold $2 in Hy, then u* itself is a Lipschitz solution to 
(1.2). Otherwise, set s = ^2 +^+(^ 2 ) choose a point (x,f) G such that 


\Du*\ > s in H X (0,t), \Du*{x,i)\ G 

Then we can do the obvious as above to obtain infinitely many Lipschitz 


solutions u to (1.2), and the proof is complete. 


8. Radial and non-radial solutions 


In this final section, we prove Theorem 1.3 on the coexistence of radial 
and non-radial Lipschitz solutions to problem (1.2) when H is a ball and uq is 


a radial function. For convenience, we focus only on Case I: Perona-Malik 
type of equations; one could equally justify the same for Case II: Hollig 
type of equations. So we assume the flux A{p) fulfills Hypothesis (PM). 

Let H be an open ball in and the initial datum uq G 

satisfy the compatibility condition 

A{Duq) • n = 0 on dVL. 
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We say that a function u defined in resp.) is radial if u{x,t) = 

u{y,t) 'ix,y E n, |x| = \y\, Vt E (0,r) (n(x) = u{y) Vx,y E |x| = \y\, 
resp.). 

We have the following. 

Theorem 8.1. Assume uq is non-constant and radial. Then there are in¬ 


finitely many radial and non-radial Lipschitz solutions to (1.2) 


Proof. The existence of infinitely many radial Lipschitz solutions to (1.2) 
follows from the authors’ recent paper 


j . We remark that these radial so¬ 
lutions are not obtained through the existence result of this paper, Theorem 

ixn 

We now prove the existence of infinitely many non-radial Lipschitz solu¬ 


tions to (1.2). In the current situation, it is easy to see that the function 
u* £ U constructed in Section 6 is radial in Qt- Our strategy is to imi¬ 
tate the procedure of the density proof in Section 7 to the function {u*,v*). 
We choose a space-time box in Qt having positive distance from the cen¬ 
tral axis {0} X [0,T] of fl-r, where is sufficiently away from A{Du*) in 
L^-sense. Then as in the density proof, we perform the surgery on {u*,v*) 
only in the box to obtain a function with membership E U 

maintained. Such a surgery breaks down the radial symmetry of tt*; hence, 
is non-radial. Note also that this cannot be a Lipschitz solution to 


( 1 . 2 ). 


Suppose there are only finitely many (possibly zero) non-radial Lipschitz 
solutions to (1.2). This forces for the non-radial function E T to be 


the L°°-limit of some sequence of radial functions in Q in the context of the 
proof of Theorem |2.2| a contradiction. Therefore, there are infinitely many 
non-radial Lipschitz solutions to (1.2). □ 
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